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Note on the Computation of the Bessel 
Function I,(x). 


The Bessel function of purely imaginary argument 
= Tm + (m + + 1) 
is of frequent occurrence in applied mathematics, statistics, and asymptotic 
number theory. The calculation of I,,(x) for m and x both large presents an 
interesting problem to the computer.! It is clear that the series for J,(x), 
(though everywhere convergent) is of little practical use for large x. Even 

the well-known asymptotic formula, designed for use with large x 


4n?— 12  (4n? — 1*)(4n? — 3?) 


+ O(e-*x-4) 
fails to be of practical use when 7 is of order x. 
Bessel functions of large integral order can, theoretically, be built up 
from the values of I(x) and J:(x) (given in tables,? or found by (2)) by 
use of the recurrence formula 


(3) = — + 


This procedure, however, affords another example in which great ac- 
curacy is needed to start with (as has been pointed out in case of J,(x) by 
Cosens [MTAC p. 99]), since more and more accuracy is lost at each step. 
The value of I99(60) = 0.0725775--- would, if computed this way, require 
about 40 decimal places in I9(60) and J,(60), and 89 applications of (3). 

NICHOLSON? has given a formula for J,(x) purporting to be “suitable for 
rapid tabulation,” when either m or x is large. The formula is not given 
explicitly (except its first term) but in terms of an infinite differential 
operator with undetermined coefficients operating on a cumbersome function 
of m and x. 

The following formula, which I have not found in the literature, is quite 
effective in case m is large and x has any real value. 


where z = (1 + x’)}, and where the y’s are the following polynomials: 


24 = 3x* — 2, 
16 = x* — 4x*, 
5760 W3(x?) = 375x® — 3654x* + 1512x? — 16, 
128 Wa(x*) = — 232x° + 288x* — 32x?, 
322560 Ws(x?) = 67599x!° — 1914210x* + 4744640x* — 1891200x* 
+ 78720x? — 256, 
192 We(x?) = 103x!2 — 4242x" + 17493x8 — 14884x° + 2580x* — 48x?, 
3440640 = 5635995x"* — 318291750x'! + 1965889800x"° 
— 2884531440x* + 1135145088x* — 99783936x* 
; + 881664x? — 2048, 
4096 Ws(x?) = 23797x'* — 1765936x"* + 15252048x"? — 34280896x"° 
+ 24059968x° — 5095936x*° + 248320«* — 1024x*. 


133 


€ 
) 
° 


134 NOTE ON THE COMPUTATION OF THE BESSEL FUNCTION [,(x) 


This formula may be derived from an asymptotic expansion of the 
solution 


y = d(In I,,(mx))/dx 
of the differential equation 
+ y/x + = + x*)/x* = (nz/x)?. 
Further polynomials y,,(#), if needed, may be found expeditiously as follows. 
Define the polynomial Q,,(#) by 
and in general by the recurrence forraula‘ 


Qn = (6m — (6m — 4)1)Qn-1 44(1 — + (1 


Then the mth term 7,, of the asymptotic series occurring in (4) may be 
written 


Tn = #(— 4n)—™ f 


In other words, if we denote the coefficient of ¢* in Qn(#) by Ax™ and define 
the polynomial ¢,,(é) by 


dn(t) = Ay™4*/(3m — 2k), 
k=0 


then 
= (— 
As a matter of fact, the polynomial 
Ym(t) = dm(1 + 


has simpler coefficients than ¢,,(#) and since ¢m(z*) = m(x?) we have given 
Tm partly in terms of x. In doing this we are following MEIssEL’ whose 
formulas for J,(nx) also involve ~m(t) for m = 6. A comparison of his 
extended results reveals four errors in his ¥¢(sec? a). For denominators 


1, 4, 192, 256, 128, 768, read 1, 4, 48, 16, 8, 48. 


This incorrect polynomial has produced an error in the corresponding term 
in the formula as quoted by Watson® 


for 3072, 768, 41280, read 192, 48, 2580. 


WaTSON has another error (not due to MEISSEL) in the polynomial where 
98720 is a misprint for 78720; noted by W. G. BickLey, Phil. Mag., s. 7, 
v. 34, 1943, p. 45. 

To illustrate the behavior of (4) for a moderate value of x we take 
x = 2 so that s = §. In this case (4) becomes: 


In I,(3n/4) = — In 3)n — — (Sx/2) — 1/150n — 99/3125n? 
— 1044097/175781250n? + 868626/48828125n* 
+ 857735827 /76904296875n5 
— 4925585088 /152587890625n° 
+ 4590654771247901 /100135803222656250n’ 
---. 
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For m = 100 these various terms are, to 20 decimals, 


15.13877 11331 89030 86048 
— 2.30258 50929 94045 68402 
— 1.03051 03088 84078 40331 
— .00006 66666 66666 66667 
-00000 31680 00000 00000 
-00000 00059 39751 82222 
-00000 00001 77894 60480 
00000 00000 01115 32887 
00000 00000 00032 28031 
-00000 00000 00000 45844 
-00000 00000 00000 01166, 


and their sum is 
11.80560 58908 83465 49--- 
Thus 
Ti00(75) = 134001.44880 18810. 
D. H. L. 


*Such Bessel functions have been used recently at the statistical laboratory of the 
University of California in preparation of certain statistical tables to appear in Annals o 
Math. Statistics. See also: J. WisHart, ‘‘A note on the distribution of the correlation ratio,” 
Biometrika, v. 24, 1932, p. 454, formula (27). 

? The most extensive tables of I,(x) are in B.A.A.S., Math. Tables, v. 6, Bessel Functions, 
part I, Cambridge, 1937, Tables VI and VIII. 

. W. NICHOLSON, ‘Phil. Mag. s. 6, v. 20, 1910, p. 938-943. 
e Q’s may be checked by the relation ‘Qn(1) = 4"; also m(m + 1) fi #”-10,.(t-*)dt 
= ed ae where B, is the kth Bernoulli number, in the notation of s 
D. F. E. MEtssgEt, Astr. Nach. v. 130, 1892, cols. 363-4. 
‘ G.N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge 1922, p. 228. 


Mathematical Tables in Phil. Mag. 


The important scientific periodical London, Edinburgh & Dublin Philo- 
sophical Magazine and Journal of Science, has been published under different 
titles since 1798, in 7 series of volumes (varying in number from 11 to 68). 
In recent years the usual print abbreviation of the title of this periodical 
has been Phil. Mag. In what follows it is proposed to list, with a few added 
notes, mathematical tables in Phil. Mag., s. 3, v. 18 to s. 7, v. 34, no. 235, 
inclusive, 1841-1943. It will be observed that more than half of the Tables 
listed were published after the last world war started in 1914, a period of 
unparalleled scientific development. Many mathematical tables have been 
already published as the result of problems arising in the prosecution of the 
present war; the number is likely to be greatly increased in the next few 
years. 


S. 3 


1. G. B. Arry “On diffraction of an annular aperture,” v. 18, 1841, p. 7. Table of 
E{= Jo(e)] = (1/2) J2* cos (e cos 0)d0, and of E*, for e = [0.0(0.2)10.0; 4D]. Airy had 
an earlier table of 2Ji(e)/e of the same scope in Cambridge Phil. So., Trans., v. 5, 1835, 
p. 291. In B.A.A.S., Math. Tables, v. 6, Bessel Functions, Cambridge, 1937, Jo is given, 
among other values, for e = [0.0(0.2)25.0; 8D]. The value of Jo(40), to 7D, was computed 
by W. R. Hamilton, Phil. Mag., s. 4, v. 14, 1857, p. 381. 
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2. J. W. Lussock, “On the determination of the numerical values of the coefficients in 
any series consisting of sines and cosines of multiples of a variable angle,” v. 33, 1848, p. 
106-127. “In the following methods of developing any function of sines and cosines by 
means of particular values of the function, the artifice by which the coefficients are suc- 
cessively eliminated is taken from Leverrier’s ‘Développements sur plusieurs points de la 
théorie des perturbations des planétes.’’’ Lubbock gives various tables for evaluating 
numerical constants when the method is practically employed. 


S. 4 


3. J. W. L. GratsHer, “On a class of definite integrals—Part II,” v. 42, 1871, p. 431-436, 
“Tables of the error function.” Values of J2°e~“"du are given for x = [3.00(0.01)4.50; 7S]. 
The range x = [0.00(0.01)3.00; 8D to 11D] was covered in a table, with first three dif- 
ferences, in C. Kramp, Analyse des Réfractions Astronomiques et Terrestres, Leipzig and 
Paris, 1799, p. 195-202. In A. A. MarKov, Table des Valeurs de l’Intégrale “dt, 
St. Petersburg, 1888, x = [0.000(0.001)3.00(0.01)4.80; 11D], and first, second, and third 
differences are given. 


S. 5 


4. H. M. Taytor, “On the relative values of the pieces in chess,” v. 1, 1876, p. 229. Table 
of “chance of the king being in check” for chessboard of n* squares (n even). 


5. J. W. Strutt, baron RAYLEIGH, “Investigations in optics, with special reference to the 
spectroscope,” v. 8, 1879, p. 266-268. For varying values of u, u~* sin* u, to 4D, is given in 


TA f= w(sinu — sin), as well as (f/fo), are given to 4D in T. II; and 


f=“ (sin w —2 sin = cos =) , as well as (f/fo)? are given in T. III. Also in Sci. Papers 
of John William Strutt, v. 1, Cambridge, 1899, p. 419-422. 


6. G. B. Arry and W. E ttis, “On a systematic interruption in the order of numerical 
values of vulgar fractions when arranged in a series of consecutive magnitudes,” v. 12, 
1881, p. 175-178. Tabular selections from their ‘Logarithms of the values of all vulgar 
fractions, with numerator and denominator not exceeding 100, arranged in order of mag- 
nitude,” Institution of Civil Engineers, Minutes of Proc., v. 65, 1881, p. 271-298; 3043 
fractions. Compare no. 7, and RMT 87. 


7. J. J. Sytvester, “On the number of fractions contained in any ‘Farey series’ of which 
the limiting number is given,” v. 15, 1883, p. 251-257. A Farey series (see Phil. Mag., 


s. 1, v. 47, 1816, p. 385) is a system of all the unequal vulgar fractions arranged in order of 


magnitude, the numerator and denominator of which do not exceed a given number, /j. 
The number of fractions in the series for the limit j, is identical with the sum of the totients 
of all the natural numbers up to j inclusive—a totient of n, denoted by r(m), meaning the 
n= 
number of numbers not exceeding m and prime to it, T(j) = z r(n). There is a table for 
n = 1(1) 500, of r(m), T(x), (3/x*)n*. In Sylvester’s Coll. Math. Papers, v. 4, Cambridge, 
1912, p. 103-109 this table is extended to m = 1000. For n = 688, for r(n) 536, read 336. This 
error is noted in J. W. L. GLAIsHER, Number-Divisor Tables (B.A.A.S., Math. Tables, v. 8), 
Cambridge, 1940, where Sylvester’s table is extended to m = 10000. See further D. H. 
LEHMER, Guide to Tables in the Theory of Numbers (Nat. Res. Council Bull. no. 105), Wash- 
ington, D. C. 1941. Compare no. 6. 


8. C. E. HoLtanp, P. R. Jones, C. G. Lams, ‘Table of zonal spherical harmonics with a 
short explanation and some illustrations of its use by John Perry.” v. 32, 1891, p. 512-523. 
Also in Phys. So. London, Proc., v. 11, 1892, p. 221-233. This contains the first table of 
P,(cos 6), m = [1(1)7; 4D], @ = 0°(1°)90°. It was reprinted in W. E. ByeRLy, An Elementary 


— _»-—-— 


7 
I 
I 
t 
j 


4 
| 
% 
| 


MATHEMATICAL TABLES IN PHIL. MAG. 


Treatise on Fourier’s Series and Spherical, Cylindrical, and Ellipsoidal Harmonics... , 
Boston, 1895, p. 278-279. A. H. H. Tatiovist (Soc. Sci. Fennicae, Acta, v. 33, no. 4, 1908) 
pointed out the numerous and non-trivial errors in this table which Byerly faithfully 
reproduced except in one case, P,(cos 7°), where he correctly put — 0.0038 for 0.0038 in 
the original. Of the 92 errors there are 38 in P «(cos 6) and 48 in P7(cos 6). Of the worst errors 
samples are as follows: P2(cos 76°), for — 0.4112, read — 0.4122; Pe(cos 7°) for 0.8476, 
read 0.8492; P.(cos 47°) for — 0.4252, read — 0.4227; P2(cos 7°) for 0.7986, read 0.8016; 
P;(cos 18°) for 0.0289, read 0.0248. A correct form of the table, including Ps(cos 8), is given 
in Tallqvist’s Grunderna af Teorin for Sferiska Funktioner . . . , Helsingfors, 1905, p. 386- 
387. 


9. J. W. Strutt, baron RAYLEIGH, “On the instability of a cylinder of viscous liquid under 
capillary force,” v. 34, 1892, p. 152-153. There is a table of x* + 1 — x*J,*(x)/I:*(x) 
= x{Je(ix)/J2(ix) + 1+ 1/x*], for x = 0.0, 0.2, 0.4, 0.6, 1.0, 2.0, 4.0, 6.0, to 3D or 4D. 
Also Scientific Papers, v. 3, 1902, p. 591-592. 


10. H. A. RowLAnD, “Gratings in theory and practice,” v. 35, 1893, p. 414. Table of “ghosts” 
(Bessel functions) J,2, for m = [0(1)14; 3D] and varying values of the parameter 0.0 to 10.0. 


11. J. Perry and H. F. Hunt, “The development of arbitrary functions,” v. 40, 1895, 
p. 508-509. Four tables of xJi(x) to 4D, for varying values of a parameter. 


12. H. NaGaoxA, “Diffraction phenomena in the focal plane of a telescope with circular 
aperture, due to a finite source of light,” v. 45, 1898, p. 9. Tables of I(/) = 1 — Je(r) 
— J,*(r) for varying values of r, to 5D. 


13. B. A. Smitu, “Tables of «Jo(x) — Yo(x) and «J:(x) — Yi(x),” v. 45, 1898, p. 122-123. 
In an article by J. H. Michell on the ‘“wave-resistance of a ship.” Tables for 
x = [0.00(.01)10.3; 4D]. «Jo(x) — Yo(x) = — sin (x sin Smith's 
tables of Yo(x) and Yi(x), for x = [0.00(.01)10.2; 4D], are given in Mess. Math., n.s., v. 
26, 1896, p. 98-101. This was the earliest table of functions of the second kind. 


14. K. PEARSON and A. LE, “On the criterion that a given system of deviations from the 
probable in the case of a correlated system of variables is such that it can be reasonably 
supposed to have arisen from random sampling,” v. 50, 1900, p. 175. Improbability 


Paew(14% for n’ odd; 
2 2-4 2-4-6---(n’ — 3) 
probability 
3 n'—3 
for n’ even. 


Table of P for nm’ = 3(1)20, x? = 1(1)10(5)30(10)70. 

S. 6 
15. W. B. Morton, “On the propagation of polyphase currents along a number of parallel 
wires,” v. 1, 1901, p. 566. = (sin x/n)—2 29"/")[sin (2x/n)]—2 © Table of » for 


gq = 1(1)6, m = 2(1)12. 
80 3 


16. H. Hitton, “A note on van der Waals’ equation” [» ak may ee v. 1, 1901, 
p. 579-589. Tables I-II for 86 = 4, 5, 6, ---16, 24, 32, 40, 48, ---80, — 8, forx = 4, ---, 
5.0; 86 = 4, — 8, for x = .3 --- .1, — .5, ---— 5.00. Graphs. There are further tables and 


graphs in Hilton's “A further note on van der Waals’ equation,” v. 2, 1901, p. 108-118. 
See also various tables in J. P. Dalton, “On the saturation constants according to van der 
Waals’ equation,” v. 13, 1907, p. 517-522. 


17. A. A. Ross, “On the conduction of electricity through gases between parallel plates 
—part II,” v. 10, 1905, p. 672-675. A table gives P, Q, R, S, T, in terms of cosh w 
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= [1.0(0.1)10; 3D]; 


S = cosh? ude; +2 


18. H. G. SavipcE, “Tables of the ber and bei and ker and kei functions, with further 
formulae for their computation,” v. 19, 1910, p. 49-58; also in Phys. So. London, Proc., 
v. 22, 1910, p. 105-114. Tables of ber x, bei x, ker x, kei x, for x = [1(1)30; 4S]; other tables, 
p. 54-57, with a similar range. SAVIDGE gave also tables of ker and kei and their first deriva- 
tives for x = [0.1(0.1)10.0; 7D-9D], in B.A.A.S., Report, 1915, p. 36-38. 


19. J. R. Arrey, “‘Tables of Neumann functions G,(x) and Y,(x),” v. 22, 1911, p. 658-663. 
Tables of Go(x), Gi(x), Yo(x), Yi(x), for x = [0.1(0.1)16; 7D]. 


20. P. F. Warp, ‘‘The transverse vibrations of a rod of varying cross-section,” v. 25, 1913, 
p. 85-106. Various tables corresponding to assumed forms of the rod; Bessel functions basic. 


21. J. R. Arrey, “Bessel and Neumann functions of equal order and argument,” v. 31, 
1916, p. 521-527. The following results are given: first root of Jiooo(z) = 0, to 2D; the 
values of J,(m) form = 6 (6D), = 48 (9D), m = 750 (14D), m = 109 (6D). Also of J,(109) 
for nm = [100(1)108; 6D]. Also of G,(m) and Gn_i(m) for n = [7(1)13; 6D]. And of G,(104) 
for m = [100(1)107; 6D]. See B.A.A.S., Report, 1916, p. 92-96 where ArrEy has elaborate 
tables of Jn-s(m), Gn(m), Gn-i(m), Yn(m), Yn-i(m) for m = [1(1)50(5)100(10)200(20)- 
400(5S0) 1000(100)2000(500)5000(1000)20000(5000)30000(10000) 50000, 100000, 500000, 
1000000; 6D]. Compare nos. 23, 37. 


22. J. R. Atrey, “The roots of Bessel and Neumann functions of high order,” v. 32, 1916, 
p. 7-14. The first five zeros of Jio(x), Ji00(x), J1000(x), Gioo(x), and Yioo(x) are given to 6S. 
ArreEy’s “The roots of the Neumann and Bessel functions,’ Phys. So. London, Proc., v. 
33, 1911, p. 219-24, has the following tables: (a) The first 40 zeros of Yo(x), Yi(x), Yo(x), 
to 5D; (b) The first 10 zeros of No(x), Ni(x), Jo(x) + Yo(x), Ji(x) + Yilx), Jo(x) — Yolx), 
to 5D. 


23. G. N. Watson, “Bessel functions of equal order and argument,”’ v. 35, 1918, p. 369. 
Incidentally there is a table of nJ/q'J,(nx)dx for n = [1(2)23; 7D]. Compare no. 21. 


24. H. A. WEBB and J. R. Arrey, “The practical importance of the confluent hypergeometric 
function,” v. 36, 1918, p. 137-141, and plate VI. Table of M(a, y, x) for y = [1(1)7; 4S], 
a =— 3,0(0.5)4.0, x = 1(1)6(2)10. 


25. J. R. Arey, “‘The addition theorem of the Bessel functions of zero and unit orders,” 
v. 36, 1918, p. 238-242. Jo(x), Go(x), Yo(x), Ji(x), Gi(x), Yi(x) are given for x = 9 and 10, 
to 14D; also auxiliary tables for the computation. 


26. J. R. Arey, “Bessel functions of small fractional order and their application to problems 
of elastic stability,” v. 41, 1921, p. 201-205. First and second zeros of J,(x) are given for 
n = [2, 3, 4, 4; 2D] (also the first three zeros for » = [— 3, — 3, — 4, — 3; 2D]), and 
n = [— 3(0.025)1; 4D], uncertain fourth decimal places, and the first zeros of J4:/3(x) 
and J41/(x), being also indicated. This last table is reprinted (without indicated un- 
certainties) in the work of Gray-Mathews-Macrobert, T. XIII, p. 317. Airey gives also the 
second zero, and (2m)~} X first zero, of J,’(x) for n = [— $(0.025)4$; 4D]; and [}(m + 1)] 
X the square of the first zero of J,(x) for » = [— 1(0.025)0; 4D]. Compare no. 27. 


27. A. Ono, “On the first root of Bessel functions of fractional order,” v. 42, 1921, p. 1021. 
With reference to Airey’s paper (no. 26), gives first zeros of J_1/2, J~3/1, J-1/3, Jo, Jis2, Ji; 
3D. 


28. H. E. H. Wrinch and D. M. Wrincu, “Tables of the Bessel functions J,(x),”" v. 45, 
1923, p. 846-849. Table I: I,(x), m = 0(1)6, x = [5(1)15; 6S]; T. II: Log T(x + 1) for 
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x = (1(1)28; 7D]; T. III: (2ex)4, e*, e~*, e=/(2xx)*, e~=/(2xx)4 for x = [1(1)20; 6S]. See also 
no. 30. 


29. F. F. P. Brsacre, “The calculation of the skin effect in electrical conductors,” v. 45, 
1923, p. 1035-1038. There are tables (independently checked) for the following functions 
in [r, 6] and [a + ib] forms: Table I: ix), ix); T. II: He(V—ix), 
Hwy — 4x); for x = [0.0(0.2)6.0; 4 or 5S]. The tables for the J’s in the [r, 6]-form were 
calculated from the corresponding tables in the (a + 1b)-form, due to W. S. Atpis, R. So. 
London, Proc., v. 66, 1900, p. 42-43. Some columns in Table II were taken from JAHNKE 
and Emp, Funktionentafeln, Leipzig, 1909, p. 139-140. 


30. H. E. H. Wrincu and D. M. Wrincu, “Tables of Bessel functions,” v. 47, 1924, p. 
62-65. Table I: I,(x), m = 0(1)6, x = [16(1)37; 6S]; T. II: (2xx)§, e*, e*/(2xx)! for 
x = [21(1)40; 6S]; T. III-IV: J2(x)/Ie2(x), and J3(x)/Is(x), for x = [1(1)15; 6D]. The 
graphs of the functions in Tables III-IV are also given. These functions are of importance 
in certain problems in acoustics; see R. So. London, Proc., v. 101A., 1922, p. 493-508. 
See also no. 28. 


31. S. P. Owen, “Table of values of the integral /o*Ko(t)dt,” v. 47, 1924, p. 736. Table for 
x = [0.02, 0.1, 0.5, 1(1)12; 6S-7S]. 


S.7 


32. H. E. H. Wrincu and D. M. Wrincu, “The roots of hypergeometric functions with a 
numerator and four denominators,” v. 1, 1926, p. 273-276. First paragraph: “During the 
last few years, the Generalized Hypergeometric Function has become of increasing im- 
portance in applied mathematics. This paper contains certain tables of the roots of certain 
generalized hypergeometric functions which have important applications in the theory of 
the vibrations of bars of various cross-sections. Attention was called to the fundamental 
significance in problems of morphology, of the periods of vibration of bars of varying 
density with cross-sections of various types, by [A.] Denby and [J. W.] NicHotson [R. So. 
London, Proc., v. 93A and 89B] in 1917.’ In this paper are given the first nine roots of the 
equation F(a; a1, a2, a3, a4; — 4(x/4)*) = 0 for a varied range of values of the parameters. 
33. H. P. MULHOLLAND and S. Go.pstEtn, “The characteristic numbers of the Mathieu 
equation with purely imaginary parameter,” v. 8, 1929, p. 839. Tables are found for the 
first eight characteristic numbers for purely imaginary g in the Mathieu equation 


dy 
a? (4a — 16gcos2x)y = 0. 


34. S. Hicucui, “On some closed algebraic curves and their application to dynamical 
problems,” v. 9, 1930, p. 186-190. Tables for solutions of the equation x* + y* = 1, for 
the area of this curve, and for the volumes of its revolution about x and y axes; etc., for dif- 
ferent values of n. 

35. N. W. McLacuvan and A. L. Meyers, “The polar form of the ker and kei functions, 
with applications to eddy current heating,” v. 18, 1934, p. 621-624. There are tables of 
No(z), ¢0(2z) (degrees), Ni(z), ¢1(z) (degrees), for z = [0.0(0.1)10.0; 5D to 8D]. 


No(z) = (ker’s + kei?z)4; = kei 
kerz’ 
Ni(z) = (ker”z + kei’z)4; = kei’ z 4 


36. J. R. Arrey, “Toroidal functions and the complete elliptic integrals,” v. 19, 1935, 
p. 177-188. The tabulation of toroidal functions which satisfy Laplace’s equation and are 
suitable for conditions over the surfaces of anchor-rings, is most conveniently and accurately 
effected, especially when the argument is large, through the complete elliptic integrals of 
the first and second kind, K and E. K(k*)=K,-In(4/h*) — Ke; E(k*) = E,-In(4/h4) + Ey 
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Table I gives K; and Kz with second differences, for h = [0.00000(0.00001)0.0001 (.0001) 
0.001(.001)0.170; 10, or more, D]; T. II gives, to 10 or more D, with second differences, 
E, and E; for h = 0.00000(0.00001)0.0001 (0.0001 )0.001 (0.001)0.100. 


37. J. R. Arrey, “Bessel functions of nearly equal order and argument,” v. 19, 1935, p. 
233-235. There are tables, to 6D, of J,(u) for w = 1(1)20; » = w — 1(0.1)u + 1. Cf. no. 21. 


dJ,(x) d #J,(x) 
236-243. Tables to 4D or 5D for x = 1(1)20; » = x — 1(0.1)x + 1. 


39. L. SILBERSTEIN, ‘On complex primes,” v. 19, 1935, p. 1104. The table gives the decom- 
Position of primes less than 1000 of the form 4” + 1 into sums of two squares, part of the 
author’s ms. table for P = 25,033. 


40. W. G. BickLey and J. NAyLer, “A short table of the functions Ki,(x), from m = 1 
cosh ugy 
to = 16,” v. 20, 1935, p. 343-347. Kin(x) = Kig(x) = Ko(x), the 
cosh" 

Bessel function of the second kind, with imaginary argument, as defined by Watson. The 
tables are for x = [0.00(0.05)0.2(0.1)2.0, 3.0; 9D]; m = 1(1)16. 

41. L. SivBerstern, “The roots of cosz = z,”’ v. 20, 1935, p. 528-531. Approximations 
to the real root, to 5D, and six complex roots, z = x + iy for x negative, are found. The 
real root, to 8D, and 12 other complex roots for x positive, were found much earlier by T. H. 
MILter, “On the numerical values of the roots of the equation cosx = x,”” Edinburgh Math. 
So., Proc., v. 9, 1891, p. 80-83. Compare Scripta Math., v. 4, 1936, p. 100. See also R. 


Cooper and J. Topp, “The large roots of cosz = az+c,”’ v. 21, 1936, p. 249-262. Compare 
no. 50. 


42. J. R. AtrEy, assisted by L. J. Comrie, “The circular and hyperbolic functions, argument 
x/v2,” v. 20, 1935, p. 721-731. There are here two tables: the first, p. 722-726, of sin (x/v2) 
and of cos (x/v2) for x = [0.0(0.1)20.0; 12D]; the second, p. 726-731, of sinh (x/v2) and 
cosh (x/v2) for the same ranges. Sin (x/v2) = v2[Ji(x) + Js(x) — Js(x) — J2(x) + Jo(x) 
+--+] and cos (x/v2) = Jo(x) — 2[Ja(x) — Js(x) + Ji2(x) — Jis(x) +---]. Compare 
Scripta Math., v. 4, 1936, p. 101-102. 


43. J. R. Arrey, “The circular sine and cosine functions, argument log.x,”’ v. 20, 1935, p. 
731-738. Bessel, Neumann, Struve, and other related functions of imaginary order iv, in 
ascending powers of the real variable x, include the factor x”, the tabulation of which 
involves cos (vlnx) and sin (vlnx). There are the following four tables: T. I: sin In(1+,) and 
cos In(1+ ,) for p = [0.000(0.001)0.020; 10D]; T. II: — sin In (1 — p) and cos In(1 — p) for 
the same range; T. III: sin Inx and cos Inx for x = 2, 5, 10, 100, 1000, x, to 13D or 14D; T. 
IV: sin Inx and cos Inx, for x = [0.1(0.1)20.0; 10D]. 


44. V. M. FaLKner, “A method of numerical solution of differential equations,”’ v. 21, 
1936, p. 624-640. In connection with the discussion of solutions of Y” — YY’ = 0, and of 
y’” — YY” = 0, tables are given for Y, Y’, Y”, Y’” for varying values of X. 


38. J. R. Arrey, “The Bessel function derivatives v. 19, 1935, p. 


45. R. C. CoLwe.t and H. C. Harpy, ‘The frequencies and nodal systems of circular 
plates,” v. 24, 1937, p. 1046. On this page are given, incidentally, two tables; in one, the 
values of the modified Bessel functions I(x), for x=0, 1, 5, 10, 20, 30; in the second the 
first ten roots of J,(x) = 0 for m = [0(1)5; 3D], 34 of 60 values incorrect; cf. MTE 21. 


o dx 
46. J. R. Arrey, “The radiation integrals f 7e-D’ 
These integrals “‘occur in a number of physical and astronomical problems, e.g., when 
n = 1, the X-ray analysis of crystal structures; when m = 2, the radiative equilibrium of a 
planetary nebula and the temperatures of stars and novae; when m = 3, temperature of 
the nuclei of planetary nebulae and the exchange of energy between monatomic gases and 
solid surfaces.” x = [0.0(0.1)10.0; 6D], a = 0, 3, 4, 2, 1, 2. Compare v. 34, 1943, p. 602-606. 


” vy. 25, 1938, p. 273-282. 
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47. A. FLetcuer, “A table of complete elliptic integrals,” v. 30, 1940, p. 516-519. “The 
following table was made in connection with work on Laplace coefficients in dynamical 
astronomy. Doubtless this is not the only case in which the modulus & is preferable to 
sink or k* as argument, and the functions of the complementary modulus &’ are not 
required.” K = x/2M. A table of K, E, and M, is given for k = [0.00(0.01)1.00; 10D]. See 
also Fletcher’s ‘““Tables of the two chief Laplace coefficients,”” R.A.S., Mo. Notices, v. 99, 
1939, p. 259-265. 


48. W. G. BickLey and J. C. P. Miter, “Numerical differentiation near the limits of a 
difference table,” v. 33, 1942, p. 1-14 + 4 folding plates. See RMT 93. 


49. C. A. Coutson and W. E. Duncanson, “Some new values for the exponential integral,” 
v. 33, 1942, p. 754-760. In T. I, Ei(x) and — Ei(— x) are given for x = [15(1)50; 10S}. 
B.A.A.S., Mathematical Tables, v. 1, London, 1931, and Matu. TaBLes Project, New York, 
Tables of Sine, Cosine, and Exponential Integrals, 2 v. 1940, give values of these functions 
for 0 < x < 15. J. P. Gram, “Undersggelser angaaende Mzngden af Primtal under en 
given Grense,” Danske Vidensk. Selskabs, Skrifter, s. 6, Naturvid. og Math., v. 2, p. 268- 
272, gives Ei(x) for x = [10(1)20; 20D], [5.0(0.2)20; 10D] with first differences. But 
— Ei(— x) in the range 15 < x < 20 had not been earlier computed. The only previous 
table of — Ei(— x) for x > 20 is that of T. Akanrra, “Tables of e~*/x and /2°e~“du/u 
from x = 20 to x = 50,” Inst. Phys. Chem. Research, Tokyo, Sci. Papers, Table no. 3, 
1929; x = [20.00(0.02)50.00; 6S, 5S] with first and second differences. Coulson and Duncan- 
son make no direct reference to the paper of Gram, from which it is clear that their last 
unit in Ei(20) should be 6, not 7. ‘‘Ei(x) is repeatedly needed in the evaluation of integrals 
concerned with nuclear structure, molecular structure and flow of heat.” 


50. A. P. HiLLMAN and H. E. Sauzer, “Roots of sin z = z,” v. 34, 1943, p. 575. In this article 
is a table of the first ten non-zero roots of sin z = 2, in the first quadrant, to 6D. Compare 
no. 41. 


Names 

ArREY 19, 21, 22, 24, 25, 26, Harpy 45 NICHOLSON 32 

27, 36, 37, 38, 42, 43, 46 Hicucui 34 Ono 27 
Airy 1, 6 HILLMAN 50 OweEN 31 
AKAHIRA 49 HILTon 16 PEARSON 14 
Apis 29 HOLLAND 8 Perry 8, 11 
BIcKLEy 40, 48 Hunt 1l RAYLEIGH 5, 9 
BISACRE 29 JAHNKE 29 Ross 17 
BYERLY 8 JonEs 8 ROWLAND 10 
CoLWELL 45 KRamp 3 SALZER 50 
ComrRIE 42 Lams 8 SAVIDGE 18 
Cooper 41 LEE 14 SILBERSTEIN 39, 41 
CouLson 49 LEHMER 7 Situ 13 
DaLTon 16 LEVERRIER 2 SYLVESTER 7 
Denpy 32 LuBBock 2 TALLQVIST 8 
DUNCANSON 49 McLacuian 35 TAYLOR 4 
ELLIs 6 Markov 3 Topp 41 
EMDE 29 MEYERs 35 WAALS, VAN DER 16 
FALKNER 44 MICHELL 13 Warp 20 
FLETCHER 47 MILLER, J. C. P. 48 Watson 23, 40 
GLAISHER 3, 7 Miter, T. H. 41 WEBB 24 


GOLDSTEIN 33 
Gram 49 
Gray 26 
HAMILTON 1 


Morton 15 
MULHOLLAND 33 
NAGAOKA 12 
NAYLER 40 


Wrincu, D. M. 28, 30, 32 
Wrincu, H. E. H. 28, 30, 32 


R.C.A. 
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121[A, C, D, Q, R]—War DEPARTMENT, Technical Manual, TM 5-236, 
Surveying Tables, July 10, 1940. Published: Washington, U. S. Govern- 
ment Office, 1940, 392 p. 13.8 X 21.5 cm. Supplement, TM 5-236, C1, 
July 29, 1942, 4 p., a list of corrections in Tables II, LIII, and LV. For 
sale by the Superintendent of Documents, Washington, D. C., $0.40. 


Of the 56 tables in this volume we shall give only brief indications of the titles of the 
first 52. These are as follows: I-II, Common logarithms of numbers, and of the circular 
functions; III-V, Natural sines, cosines, tangents, cotangents; VI, Stadia reductions; VII, 
Elevations above sea level from barometric readings; VIII, Coefficients for temperature 
corrections, barometric; IX—X, Conversion of meters to feet and feet to meters; XI, Con- 
version of seconds of latitude and longitude to feet; XII, Tables of S; XIII, Conversion of 
decimals of a day to hours; XIV, Inclination corrections for 50 meter tape; XV, Curvature 
and refraction; XVI, Corrections for vertical angles; XVII, Common logarithms of m, 
spherical excess factor; XVIII, Azimuths of Polaris at elongation, 1940-1950; XIX, Con- 
version of arc into time, and time into arc; XX, Standard time signals; XXI, Mean astro- 
nomic refraction, with corrections; XXII, Parallax and altitude of the sun; XXIII-XXIV, 
Conversion of mean time to sidereal time, and of sidereal time to mean time; X XV, Arc-sin 
cerrections for inverse position computations; XX VI, Conversion of degrees Centigrade to 
degrees Fahrenheit, and degrees Fahrenheit to degrees Centigrade; XX VII, Common loga- 
rithms of versed sines and external secants; XXVIII, Natural versed sines and external 
secants; XXIX, Methods of expressing gradients; XXX, Differences in elevations for given 
horizontal distances and gradients from 0° to 30°; XXXI-XXXII, Differences in eleva- 
tions, and horizontal distances, for given slope distances and gradients from 0° to 30°; 
XXXIII, Reduction to the horizontal of distances paced on slopes; XXXIV, Radii and their 
logarithms; XXXV, Tangent offsets and middle ordinates; XXXVI, Tangents and ex- 
ternals for 1° curves; XX XVII, Corrections for tangents and externals; XX XVIII, Long 
chords and actual arcs; XX XIX, Corrections for subchord lengths; XL, Circumferences 
and areas of circles; XLI, Squares, cubes, square roots, cube roots; XLII, Conversion of 
inches and sixteenths to decimals of a foot; XLIII, Conversions of sixteenths of an inch to 
decimals of an inch; XLIV, Map scales in English measure and metric units; XLV, Minutes 
into decimals of a degree; XLVI, Useful constants and formulae; XLVII, The Greek al- 
phabet; XLVIII, Lengths of circular arcs for unit radius; XLIX, Corrections to AY for 
magnification of scale; L, Reduction of geographic azimuth to grid azimuth; LI, Ageton 
table for meridian determinations; LII, Common logarithms of circular functions in 
hundredths of a degree. 

The next four tables involve mils. Table LIII gives the logarithms, to 5D, of sine, 
cosine, tangent, cotangent, for every mil of the quadrant. Tables LIV-LV give natural 
values, to 5D, of these same functions, for the same range. Table LVI is for conversion of 
degrees to mils, and mils to degrees. 

The supplement gives one correction in Table II, and 48 in Table LIII. Also two pages 
to replace two of the four pages of the table of Natural tangents and cotangents in mils. 

The mils table calculated by ANDOYER to which reference has earlier been made (MTAC, 
p. 39, 84) may be noted once more, Service Géographique de |’Armée. Tables de Logarithmes 
a cing Décimales pour les nombres de 1 & 12 000, et pour les Lignes Trigonométriques dans le 
Systéme de la Division de la Circonférence en 64000 parties égales (dixiémes du milli¢me de 
lV'artillerie). Paris, Imprimerie Nationale, 1916; with differences and proportional parts. 

See also the description (MTAC, p. 40) of punched cards for the six natural trigonometric 
functions, to 5D, for every tenth of a mil in the quadrant. Also RMT 126. In 1942 A. W. 
TuckKER published a sheet of four-place tables, for each 10 mils, of the six logarithmic and 
natural trigonometric functions. 
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122[A, D].—War DEPARTMENT, Technical Manual, TM 4-238, Coordinate 
Conversion Tables March 25, 1943. Printed at the U. S. Government 
Printing Office, Washington, D. C., 1942. 338 p. 13.8 X 21.5 cm. The 
tables in this volume were computed by the MATHEMATICAL TABLES 
Project, New York. For sale by the Superintendent of Documents, 
Washington, D. C., $0.40. 


This volume contains merely a single page of Introduction; Table I, coordinate con- 
version tables (p. 2-321); and Table II, eight-place table of sines and cosines of angles 
expressed in mils (p. 322-337). Table I gives the values of A sinx and A cosx, where x is 
in mils, a mil being defined as one-sixteen hundredth of a right angle. A = 1(1)100 units, 
50 units (and 10 mils) to a page. With left-hand pages giving Asinx and right-hand pages 
Acosx, four pages are therefore required for each 10 mils, and hence 320 for the quadrant. 
The table is to 4D. If the straight line to a point 7654 yards distant make an angle of 81 
mils with a line of reference a coordinate of the point, such as Acosx, is found as 7629.81 
(7600cos81 + 54cos81 = 7575.98 + 53.83). The result will be correct to within a unit in 
the second decimal place. When higher accuracy is required, Table II may be used, since 
this gives, to 8D, the values of sines and cosines, with first differences, for every mil in the 
quadrant. Identically the same second differences are given at the bottom of each of the 
pages 322-337 of the table, whereas it is only correct for page 332. This fact is undoubtedly 
due to the unscientific methods of the Government printer rather than to the personnel 
of the MATHEMATICAL TABLES Project. Note added in proof: The necessary correction 
sheets (TM 4-238 C1, 3 p.) were issued in August, 1943. 


123[A, F, P].— Formulas in Gearing with Practical Suggestions. Seventeenth 
ed., Second printing, Sept. 1943, Providence, R. I., Brown & Sharpe 
Mfg. Co., 15.1 X 22.9 cm. $1.50. 


Since details of this excellent work have already been given (see RMT 88), it will now 
simply be noted that in this edition a correction has been made on p. 224. The gear ratio 
of the decimal equivalent .4737 is 9/19, instead of 9/10 given in previous editions. 


124{A, N].—A. A. Nikitin, Tabliisy Vychislenita Protsentov (Tables for cal- 
culating percentages], second ed., edited by P. P. Andreev, Moscow, 
1939, 733 p. 21.5 X 13.9 cm. 18.90 roubles. 


The prefaces to the first and second editions, by Andreev, occupy p. 3-7, while a state- 
ment of the purpose and description of the tables is on p. 8-13. The Table (p. 14-733) is 
for percentages 1%(1%)240%, and for numbers 100(1)999. Results are given to 1D. Each 
page of the Table is arranged in 30 horizontal lines (for percentages) and 10 vertical columns 
(for numbers), so that the whole range of 240 percentages occupies 8 pages for each ten 
numbers, and therefore the complete set of tables occupies 90 X 8 = 720 pages. Moreover, 
on each page there are two tables of proportional parts, for percentages on the left margin, 
and for numbers on the right margin, the latter being subdivided into three sections, for 
10 horizontal lines each. There were 10 000 copies in this edition. 

The volume may obviously also be viewed as another Russian table for obtaining the 
product of two numbers. Compare N. 16. 

S. A. Jorre 


125(C, D].—J. T. Peters, Siebenstellige Logarithmentafel. 1: Logarithmen 
der Zahlen, Antilogarithmen, Additions- und Subtractionslogarithmen nebst 
einen Anhang mit Formeln und Konstanten; and II: Logarithmen der 
trigonometrischen Funktionen fiir jede zehnte Sekunde des Neugrades, log 
sin und log tg von 0*.0000 bis 3*.0000, sowie log cos und log ctg von 97.0000 
bis 100#.0000 fiir jede Sekunde (1°° = 0*.0001) des Neugrades. Berlin, 
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Reichsamt fiir Landesaufnahme, 1940, 2 v.; vii, 493, and vi, 666 p. 
17.2 X 26 cm. 33.50 marks and 45.00 marks. 


The second volume of mathematical tables by Peters was (in collaboration with Bau- 
schinger) an Eight-place Table of the Logarithms of all Numbers from 1 to 200 000 (1910, 
revised ed., 1936); proportional parts were given throughout. The three principal eight- 
place tables of logarithms of numbers previously published were those of (1) JoHN NEWTON, 
London, 1658, for numbers up to 100000; (2) SERvicE G&oGRAPHIQUE DE L’ARMEE, 
abridged from the Tables du Cadastre (see MTAC, p. 34, 85), Paris, 1891, up to 120 000; 
(3) MenpizABaL TAMBORREL, Paris, 1891, up to 125000 (see MTAC, p. 40). The Bau- 
schinger-Peters volume naturally superseded all of these. This table was followed in 1922 
by a work of Peters alone, namely a Ten-place Table of the Logarithms of all Numbers 
from 1 to 100 000. From this last volume the present Seven-Place Table has been prepared. 

Except for 99 numbers the abridgement was readily effected. For 96 of these numbers 
the correct seven-place values were determined by consulting a twelve-place manuscript 
table containing the logarithms of numbers from 1 to 100 000 with an error of at most a 
unit in the twelfth place; this manuscript was already used in preparing the above men- 
tioned work of 1910. The three remaining doubtful cases were recalculated and a table of 
the logarithms of the numbers 1 to 100 000 was compieted, with no error in the seventh 
place as much as 0.5. The formula used by Peters for calculating the logarithm of primes, 
was the rapidly converging series 


1 1 

While this table will be of value it should be borne in mind that no seven-place table 
has been published to supersede the work of 73 years ago by EDwARD SANG, A New Table 
of Seven-Place Logarithms of all Numbers from 20000 to 200 000, London, Layton, 1871 
(final printing 1915). Its preparation was based on a Sang fifteen-place manuscript of the 
logarithms of all numbers to 370 000, now in the Library of the Royal Society of Edinburgh. 

The second section (p. 187-387) of the first volume under review contains a Seven-place 
Table of the Antilogarithms of numbers or seven-place numbers for the mantissas 00 000 
to 99 999. As the basis of his calculations Peters used the work of HOLGER Prytz, Tables 
d’ Anti-Logarithmes, published by the Royal Academy of Sciences of Copenhagen in 1886. 
It contains fifteen-place antilogarithms for mantissas 000(001)999. Peters finally developed 
a twelve-place antilogarithm table, with mantissas 00 000(00 001)99 999, in which the 
error in the seventh place was always less than 0.502. In preparing his seven-place table 
with an error in the seventh place less than 0.5 in every case, comparison was made with 
R. Shortrede’s seven-place table of antilogarithms which is the second table in his Logarith- 
mic Tables to Seven Places of Decimals, Edinburgh, 1844. There were other editions of this 
table but none of them, apparently, in the year 1854 which Peters gives. 

The prefaces of Peters often contain interesting material, and the present one is no 
exception when he introduces matters now presented. One of the principally used logarithmic 
tables making possible the determination of both the logarithm of a number and the aumber 
corresponding to a given logarithm, exhibits also a table of antilogarithms in the solution of 
both problems. Why, it may be asked, does one then still calculate a table of antilogarithms? 
As principal reason for this, Peters remarks, we present the following: The differences of a 
table of antilogarithms are only about half as large as the differences of an ordinary table 
of logarithms. H. Bruns writes about this in his Grundlinien des wissenschaftlichen Rechnens 
(Leipzig, 1903) p. 49: “For these reasons it were indeed very much more to the point to 
tabulate, not the logarithm but the inverse function ‘Number x = 107.’ In the usual seven- 
place arrangement for log x the first differences run from 43 to 435, and the group of three- 
place differences amounts to 48%. In the corresponding tabulation of the Number one would 
have on the contrary the differences running from 22 to 231, and the total of the three- 
figure differences is less than 37%, while the extent of the table has only increased from 180 
to 200 pages.” A further advantage of employing tables of antilogarithms lies in the fact 
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that one may more readily determine the number corresponding to a given logarithm with 
an antilogarithmic table than with a logarithmic. 

The third section (p. 389-474) of the first volume is devoted to Tables of Addition and 
Subtraction Logarithms. For the calculation of these the fundamental tables were those of 
M. H. Andoyer, ‘“Tables fondamentales pour les logarithmes d’addition et de soustraction,”” 
Bull. Astronomique, s. 2, v. 2, 1922, p. 5-32. They give sixteen-place values of A and S, 
for D = 0.00(0.01)9.00. (Compare RMT 116). With this as a basis Peters developed ten- 
place tables for A with D = 0,0000(0.0001)6.94; and for S with D = 0.3000(0.0001)6.94. 
The corresponding seven-place tables made from these had an error less than 0.502 in the 
seventh place. The final checking to assure that no error in the seventh place was as much 
as half a unit was carried through by G. Witt, who had made an important contribution 
to the Anhang of Peters and Stein in Peters’ Zehnstellige Logarithmentafel, v. 1. 

If D = loga — log}, and loga > logb, from the A table, p. 389-434, is found 
log (a + 6) = loga + A. 

In connection with the S table, p. 435-474, there are two cases. When D > 0.3000, 
log (a — 6) = loga — S. If, however, log a — log b < 0.3000, set log a — log b = S and 
find the D corresponding to S in the S table, then is log (ac — 6) = loga — D. 

At the top of every page of the first of the four tables to which we have referred is log N; 
of the corresponding pages of the second table, N; of the third table, A; and of the fourth 
table, S. In the right-hand columns of all of the pages are differences, d, and at the bottoms 
of the pages, proportional parts. 

The final section of the first volume (p. 475-493) is devoted to an admirable collection 
of Formulas and Constants [goniometry (higher trigonometry); series; plane and spherical 
trigonometry; solutions of quadratic, cubic and linear simultaneous equations; theory of 
errors]. 

The second volume contains seven-place tables of the logarithms of the trigonometric 
functions, sine, tangent, cotangent, cosine, for each 0*.001. (Compare MTAC, p. 33-39). 
As basis for preparation of this table Peters used his Ten-place Table of the Logarithm 
of the Trigonometric functions from 0° to 90° for each 0°.001. Since 0*.01 = 0°.069 he 
took every ninth value from the ten-place table and then inserted in each interval 0*.01 
nine new values at intervals 0*.001. Thus a ten-place table was made for values of log cos 
0*.000 to log cos 50*.000, log tan 3*.000 to log tan 50*.000 and log sin 3*.000 to log sin 50®.000 
for each thousandth of a grade, or for each tenth centesimal second. The range 0*.000 to 
3.000 receives special treatment where tables of the usual S = log (sinw/w) and 
T = log(tanw/w) are introduced for determining log sinw = S + logw and log tanw = T 
+ logw, w being in grade units. The preparation of the ten-place table of S and T called 
for twelve-place calculation. The still missing value for each log cot was at once written 
down by simply subtracting the corresponding value of log tan from ten. 

The next step was to derive from this manuscript a seven-place table in which the 
error in the seventh decimal would be at most 0.005 of a unit. In seeking out all the doubtful 
values, comparison was made with the eight-place table derived from the tables du cadastre, 
namely: Tables des Logarithmes @ huit Décimales des Nombres entiers de 1 2 120 000 et des 
Sinus et Tangentes de dix Secondes en dix Secondes d’Arc dans le Systeme de la Division 
Centésimale du Quadrant, Paris, 1891. In this extraordinarily correct table Peters found only 
two errors, the first of which he states as follows: log sin 25*.830 = 9.596 26245, should be 
9.596 26245", so that log sin 25%.830 = 9.596 2624. The second error noted, in log cot 34®.536, 
is one of those published in 1891 by Mendiz4bal Tamborrel, and reprinted in MTAC, p. 85. 

This table occupies p. 136-665. Proportional parts are at the bottoms of the pages 
throughout. The pages 1-133 are occupied with a seven-place table of logsin w and logtan w, 
for the range w = 0*.0000 to w = 3.0000, for each centesimal second. G. Witt was again 
a calculator for this volume. It does not seem as if the typographical display of either 
volume could have been improved. 

The preface of the first volume is dated Berlin, February, 1940; and of the second, July, 
1940. These volumes were received at the Copyright Office of Library of Congress on June 
12, and December 3, 1940, respectively. There are film copies in the Library of Brown 
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University. We have yet to learn whether or not the marvellous Peters published one or 
more other volumes in the final year of his life. 
R.C. A. 


126{[D].—AprLieD MATHEMATICS PANEL, National Defense Research 
Committee, Tables of Trigonometric Functions to five significant figures 
and for every tenth of a mil, AMP Report no. 24.1, September, 1943. 
x, 320 p. 21.2 X 27.4. Printed by the photo-offset process from manu- 
script. This volume is only available to members of the Panel. 


This table for every tenth of a mil in the quadrant (6400 mils = 360°), gives the values 
for sine (5D), tangent (5D), cotangent (4D), and cosine (5D), on left-hand pages, and for 
secant (5D), 4 tangent (5D), 4 cotangent (4D), cosecant (4D) on right-hand pages. The 
table is based on original calculations and is the most extensive table yet published of the 
natural trigonometric functions with the mil as argument; compare RMT 121, 122. Such 
tables are of importance in ballistics. 

The preparation of this table was started with values of trigonometric functions, 
taken from Peters’ Seven-Place Values of Trigonometric Functions for every Thousandth of a 
Degree (see RMT 79) at intervals of 4 mils = .225 degrees from 0° to 45°. With interpolations 
a five-place table was obtained for sin x, cos x, tan x, and cot x, at intervals of one mil. 
The values of tenths of mils were obtained by linear interpolation from the values at 
intervals of one mil. Among formulae used was secx = cosx(1 + tan*x). In making a final 
revision the Mathematical Tables Project of New York computed sin x and cosx for 
x = [0(1)200; 8D] where x is expressed in mils. These were then subtabulated (using linear 
interpolation) for tenth-mil intervals. Values of cot x, } cot x, and csc x were then obtained 
directly by division. These five functions were then rounded off to five decimal places and 
punched into International Business Machines cards. And so on, with numerous checks, 
I. B. M. equipment being used throughout. 


127[D].—Hans Hor, Seven Place Full Natural Trigonometric Tables, [v. 1} 
Sine, [v. 2] Tangent, [v. 3] Cosine, [v. 4] Cotangent, Ann Arbor, Mich., 
Box 470, University Research Assoc., 1943. Each volume contains 184 p. 
(including the title page). 21.1 X 27.9 cm. $4.00 per v. The v. numbers 
indicate merely the order of publication; v. 4 will not be published be- 
fore next month. 


On the back of each title page is “copyright 1943 By Hans Hof,’’ and on the opposite 
page of each volume is the following: ‘‘All values have been checked carefully and found 
correct within approximately one point in the last or seventh place. To insure this and the 
many values given, this book has been lithoprinted from the original. The author and 
publisher will appreciate any notification of errors.’’ In no volume is there any other text. 

The three volumes contain 552 pages of tables of the natural sines, cosines, and 
tangents, to 7D, for every second of the quadrant. Thus each page is devoted to 15 
minutes, and the table runs from the bottom to the top of each page, the volume being 
turned sideways for use. Under each minute are the values of the functions for 0” followed 
by twelve solid groups of 35 figures for each 5 seconds. Exceptions to this statement occur 
in the first five or six degrees of the sine and tangent volumes where only some of the 
required zeros were typed, or inserted by hand, in the manuscript reproduced; there is no 
uniformity of procedure in this regard. For the most part, in the sine and cosine volumes, 
the decimal point occurs only before the first and last numbers of each column, but where 
zeros are inserted at random, a decimal point is usually placed before them. In the tangent 
volume decimal points are scattered indiscriminately over the page. There are a number of 
cases in top and bottom lines where the decimal points are missing. There are no rulings and 
the top and bottom margins of each page are less than half an inch wide. In general the 
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printing is clear but there is great variation between dark and light on different pages, or 
even on the same page. Not a few figures are blurred or indistinct, but the context usually 
indicates what they must be. One of several places where this is not the case (in the reviewer's 
copy at least) is at sin 66° 24’ 11”. 

The arrangement of the sine-cosine volumes, for example, is the very inconvenient one 
of Gifford’s table, even though this form is especially in demand in optical establishments. 
Throughout these two volumes, the sines always occur at the tops of the pages, 0°-89°, 
and the seconds, 0’’-60” in the extreme left-hand columns. Thus the sines and cosines of 
the same angles must be sought in different places. The curiously titled v. 1 and v. 3 will 
naturaily be bound together. 

The author has not responded to the editors’ request for information as to his method 
of computing or compiling his tables. Do they involve any original calculations whatever? 
Or were they, perhaps, the result of merely rounding off values in volumes of GrrrorD (see 
MTAC, p. 9), or in the volume of PETERS (see RMT 78, and 128)? There is evidence against 
the Peters surmise, since a test sample led to so many differently rounded values. 

This is the first printed seven-place table of these functions for every sexagesimal second 
of arc. During the last five years at least, the British Nautical Almanac Office has been in 
possession of a completed manuscript of such tables, with the cotangent, prepared by 
PETERS and Comrie. It is most earnestly hoped that the Director of the Nautical Almanac 
may see his way clear to arrange for publication of this manuscript at an early date. 

Hof’s work was published in this country almost simultaneously with the American 
edition of the Peters-Comrie volume of Eight-place Table of the Trigonometric Functions for 
every sexagesimal Second of the Quadrant (see N 6 and RMT 128). The latter publication is 
of proved accuracy to the last decimal place, and it is admirably arranged and printed. 
Unit errors in the seventh places of Hof’s tables are readily shown to be extremely numerous, 
the whole work is very poorly printed and ill arranged, and the price is decidedly excessive. 

R.C. A. 


128[D].—J. T. Peters, Eight-Place Table of Trigonometric Functions for 
every sexagesimal Second of the Quadrant. Achtstellige Tafel der trigono- 
metrischen Funktionen fiir jede Sexagesimalsekunde des Quadranten. . . . 
Published and distributed in the Public Interest by Authority of the Alien 
Property Custodian under License No. A-152. Berlin, Landesaufnahme, 
1939. Ann Arbor, Michigan, Edwards Brothers, 1943. xii, 901 p. 21.2 
X 27.9. Photo-lithoprint reproduction. $20.00. 


The original of this work, as well as the British War Office reprint, has been already 
reviewed in MTAC, p. 11-12, 65. At the latter reference, last April, we urged “most strongly” 
that the Alien Property Custodian arrange for the reproduction of this great work of Peters 
and Comrie. For, it would then be available not only for scientific workers in this country 
but also for our British friends, since copies of their War Office reproduction can not be 
procured by those not in Government employ. It is therefore with satisfaction that we greet 
this October publication. 

The 901 pages of the table are reproduced exactly as in the original. The new title-page 
is indicated above; to the copyright notice on the back of the title-page has now been added, 
“Copyright vested in the Alien Property Custodian, 1943, pursuant to law.” Instead of the 
original German we now have an English Preface and an English Introduction, each 
occupying a single page. These exact translations from the German were prepared by Mr. 
L. J. Comrie, although there is no reference to this fact in the volume; the wording is 
naturally quite different from that of the corresponding parts of the War Office edition. 

The type page and side margins are just as large as in the original, but the height of 
the volume has been reduced by more than half an inch. In turning over the 900 pages the 
reviewer did not notice more than about a score of places where the figures were blurred 
or indistinct, and in not more than four of these cases was it impossible to guess what the 
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numbers must be. The green buckram binding, with inked title on the back, is stronger 
than that of the original. There were 1000 copies in the present edition. 

Of course the publishers do not intentionally mislead when they advertise that the 
work was published in 1939 at $24.00, because it was priced at 60.00 Marks. Any member 
of the American Mathematical Society, however, and many others, could have bought this 
volume for $15.00, plus transportation charges. Of the four copies of this work which I 
know of reaching America, the cost, apart from the one sent to Library of Congress for 
copyright purposes, was less than $18.00 a copy. The cost of reproducing such a large work 
must be considerable. But might not the return to the publishers have been greater if it 
had been listed at $12.00? However that may be, there can be no doubt as to the importance 
of having this publication available at the present time. 


R.CA., 


129{I, K].—R. L. ANDERSON and E. E. Houseman, Tables of Orthogonal 
Polynomial Values Extended to N = 104, Ames, Iowa, Iowa State College 
of Agriculture and Mechanic Arts, 1942. 14.4 X 22.2 cm. (Research 
Bulletin 297, p. 593-672, Agricultural Experiment Station, Statistical 
Section.) 


The orthogonal polynomials considered in this Bulletin are 


1 


3n? — 7 
= Note = at - = = = a}, 


(1) 12 0 

S(n? — 7 15n* — 230n? + 407 


where x is a variable which takes the values 1, 2, 3, ---, , and the \,, which depend on n, 
(except for A» = 1), are so chosen that the values of the é,’, (r = 1, 2, ---, 5), corresponding 
to the 2 values of x are integers reduced to lowest terms. For example, taking » = 4, the 
values of £3 corresponding to x = 1, 2, 3, 4 are — 3/10, 9/10, — 9/10, and 3/10, respec- 
tively, so that by choosing \3 = 10/3 the corresponding values of £3’ are — 1, 3, — 3, and 
1, respectively. It is clear that A; = 1 or 2 according as n is odd or even. If x takes any se- 
quence of equispaced values x1, X2, Xs, --*, Xn, then & can be defined by 


and the é,’ are then given by equations (1) in terms of this variable. It is evident that x 
can be assumed to take the values 1, 2, ---, m without loss of generality. 

The principal use of these orthogonal polynomials is in fitting a polynomial in x by the 
method of least squares. Let 1, yo, ---, Yn denote the values of a variable y corresponding 
to the values 1, 2, ---, , respectively, of a variable x, and let these values of y be of equal 
precision, that is, let « be the standard deviation of y, for x = 1, 2, ---, m. If the polynomial 


(2) Y = ao + aix + +---+ (k <n — 1) 


is to be fitted to these data by the method of least squares, the a’s can be evaluated by 
solving simultaneously the k + 1 normal equations 


(3) xiy, = do + a; ay (Gj = 0,1, 2, R) 
=z z z 


where 2, denotes summation over x = 1, 2, ---, 2. The sums of powers up to the 2k*® of 
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the natural numbers are required, and the solution of these equations can be very laborious, 
especially when k > 3 and n is large. However, if Y is expressed as 


(4) Y = Aoto’ + +---+ Arts’ 


where £0’, £1’, &’, ---, &&’ are polynomials of degree 0, 1, 2, ---, & in x, respectively, and 
possessing the orthogonal property 

=0 for cts 
(5) 

(é-’)? + 0, 


then the normal equations for the A’s, namely, 
(6) z Ao z + Ai Ax (r = 0, 1, 2, k), 


reduce to 
(7) z (r 0, 1, 2, k), 


and the determination of the A’s is greatly simplified, especially when the values of the ¢’’s 
are known for x = 1, 2, ---, m: The polynomials £9’, £1’, ---, &s’ defined in (1) are known 
to possess the orthogonal property (5), and the Bulletin under review makes available 
more extensive tables of their values, and of the values of 2, (,’)*, than have been avail- 
able heretofore 

The first fifteen pages of this Bulletin are devoted to methods of calculation and use, 
and to references cited. The tables begin on p. 610, and for m = 3(1)104 give the exact 


2 1 
values of £1’ to 5’ (or, to £,’ for m = 5) corresponding to x = we (1)n, or x = wo (1)n, 
according as m is even or odd; when  < 9, the full set of values for x = 1(1)n is given. 
1 
This abridged tabulation is made possible by the fact that for x = — + m the values 


of 2 are equal, and also of &; similarly &,’, &s’, and &’ assume values which are of opposite 
sign. The value of 2, (£,’)? and of A, is given for each value of m and r below the corre- 
sponding column of values of £,’. 

The portion of these tables for » = 3(1)52 is a reproduction, with some variation of 
format, of Table XXIII in R. A. Fisher and F. Yates, Statistical Tables for Biological, 
Agricultural and Medical Research, Edinburgh, 1938, and contains the erratum noted in 
MTE 9, p. 86, namely, the sum of (£2’)* for m = 39 should read “4,496,388” and not 
“496,388” as tabulated. This error has been corrected in the second edition (1943) of the 
source table. 

P. L. CHEBysHEv! considered the determination of such polynomials for the general 
case where the y’s are of varying weight and x takes any set of m distinct values. For the 
case where x takes any m equispaced values Chebyshev listed* polynomials go, ¢1, ¢2, ¢2, 
¢, ¢5, satisfying (5) and gave? a recurrence formula relating g, to gr_1, ¢r—2, and g; for 
r > 2. In consequence polynomials satisfying (5) are sometimes referred to as Chebyshev 
polynomials, but this terminology is reserved, by specialists in orthogonal polynomials, for 
another set of orthogonal polynomials which he investigated also.* 

Much later the subject was treated independently by F. Esscuer,‘ by C. Jorpan,® 
and by R. A. FisHER.® FISHER gave explicit expressions*’ for polynomials , £1, ---, & 
when the y’s are of equal weight and x takes any equispaced values, which become the 
expressions in braces in (1) when x takes the values 1, 2, ---, . The coefficients of x, in £, 
are always unity. Miss F. E. ALLAN® pointed out that FisHER’s é, are related to CHEBYSHEV'S 
¢r by 


r! 


(8) = Qn)! er 


(r = 0,1, 2, 
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Making this substitution in CHEBYSHEV’s recurrence formula? and advancing the subscript 
to r + 1 it follows that 
— 

(9) = fife — 447 — 1) 1) | 
Miss ALLAN,?® using equation (J1) below, has provided a simple proof of this recurrence 
formula. It om been used in the preparation of this Bulletin to find the values of £4: for 
x = 1, 2, ---, from the corresponding values of and The value of was then 
determined and the values obtained from £'r41 = 

Starting from CHEBYSHEV’s expression® of y, as proportional to the r*® difference of 
the product 


x(x — 1)(x — 2) --- — n)(x —n—1)--- (x +1), 


Miss ALLAN derived two general formulae (p. 316 and p. 319) for &, the latter being given 
in the third and subsequent editions of FisHER’s book,’ and she tabulated the polynomials 
, up to £10. The corresponding general formulae have been derived by A. C. AITKEN” for 
the polynomials 


(10) & = ¢r/(r!)*, (r = 0, 1, 2, -- *), 
the values of which are always integers for x = 1, 2, 3, ---, m. AITKEN also gave procedures 
for building up the polynomial values from terminal values and differences. These pro- 
cedures have been used by the authors ef this Bulletin to check the values of ¢, obtained by 
means of the recurrence relation (9). 

CHEBYSHEV’ obtained a formula for 2, ¢,*, which in terms of the £, becomes 


(r!)4 

(2r) + 1)! 
(r!)* 

~ (2r) + 1)! 


This formula is derived directly by Miss ALLAN,’ and the corresponding formula for the T;, 
is derived by AITKEN.” A modification of this formula was used by the authors of this 
Bulletin to compute the values of 2, (£,’)? given in the tables. 

The reviewer has been surprised to notice that the authors of this Bulletin do not even 
mention CHEBYSHEV, whose contributions to the subject are clearly indicated in the papers 
by AITKEN and ALLAN which they cite. 


(11) (n+ 1)(n +7 —1)--- (w+ 1)n(m — 1) --- (n — 1) 


n(n? — 1?)(n? — 2?) --- (n? — 1°), (r = 0, 1, 2, ---). 


C. EIsENHART 


1 P. L. CHEBYSHEV, “‘Sur les fractions continues,” Jn. d. Math., s. 2, v. 3, 1858, p. 289- 
323. Oeuvres, St. Petersburg, v. 1, 1899, p. 201-230. 
‘Sur l’interpolation par la méthode des moindres carrés,”” Akad. Nauk, 
émoires, s. 7, v. 1, 1859, p. 1-24. Oeuvres, v. 1, p. 471-498 
gn Orthogonal Pol omials, Am. Math. So., Coll. Publ., v. 2, New York, 1939. 
iF ESScHER, ‘ ‘Ueber die — in Schweden, 1866-1914, a Lund, Observatoriet, 
Meddel., s. 2, v. 23, 1920, p. 10- 
5 CE. Jorpay, ‘Sur une ‘de polynomes dont chaque somme partielle répresente la 
meilleure approximation d’un degré donné pent la méthode des moindres carrés,”’ London 
Math. So., Proc., s. 2, v. 20, 1920, p. 297-325 
®R. A. FISHER, “Studies i ws Cro Variation.—I. An examination of the yield of dressed 
grain from Broadbalk,” Jn. A 4., v. 11, 1921, p. 107-135. 
7R. A. FisHer, Statisti "Methods ‘or Research Workers, Edinburgh, third ed., 1930. 
8F, E, ALLAN, “The general form of the orthogonal polynomials for simple series, with 
proofs of their simple properties,’ ’ R. So. Edinburgh, Proc., v. 50, 1930, p. 310-320. 
“Sur Vinterpolation,” Akad. auk, Zapiski, v. 4, 1864. Oeuvres, v. 1, 
p. 
10 A. C. AITKEN, “On the graduation of data by the orthogonal polynomials of least 
squares,” R. So. Edinburgh, Proc., v. 53, 1933, p. 54-78. 


| 
4 
7 
q 
z 
} 
| 
| 
| 
| 


<9 8 


et, 


RECENT MATHEMATICAL TABLES 151 


130[I, K]—Truman (1884- ), The Kelley Statistical 
Tables, New York, Macmillan, 1938, v, 136 p. 21.5 X 27.8 cm. $4.50. 
This book contains six different statistical tables, of which the first is the largest and 

most important. 

Table I. The area (p) under the normal curve is given by 


(1) 
and its ordinate (z) by 
(2) 


Equation (1) may be inverted to yield x as a function of p. This inverse function x(p) is 
called by statisticians the normal deviate corresponding to the probability level p. If x(p) 
is substituted in equation (2) the ordinate z becomes a function, 2(p), of the area p to which 
it is the boundary ordinate. 

The table gives, in parallel columns with argument #, the functions x(p), z() and the 
algebraical expressions vba, Vi — pand V1 — ¢, where g = 1 — p. All five functions are 
tabulated for p = [0-5000(0-0001)0-9999; 8D]. The complementary argument g (which 
statisticians call the complementary probability level) is given in the right-hand column. 
No differences are given. The table covers 100 well-printed pages of convenient size; the 
layout is good, although we should have preferred the eight decimals divided into the 
more conventional groups of five and three, rather than four and four. 

Most tables of the normal curve—e.g. (a), (c), (h) and (1)—give p and z as functions of 
x. The inverted form of normal deviates x(p) is fairly common in statistical usage; see, 
for example, (b), (g) and (j), but the ordinate z as a function of the probability p has little 
or no application, although it is given in (j). Konpo and ELDERTON’s 10-decimal normal 
deviates x(p), which are also given in (j) for p = 0-500(0-001)0-999, formed the starting- 
point of the present tables. Both (j) and KELLEy’s Tables go beyond the present-day re- 
quirements of users; four-decimal accuracy in x and 2 is all that is wanted for most applica- 
tions. The fine interval of 0-0001 in is certainly a convenience for some problems such as 
the calculation of probits (d)! in toxicology, but for statistical routine work of the nature 
of tests of significance a three-decimal (if not a two-decimal) argument is all that is required; 
the one-page tables given in (b) and (g) are handier to use. 

The functions Vq , Vi — # and V1 — ¢ are entirely independent of the columns x 
and z. The author states that he has included them in the same table simply because the 
requisite arguments were available. The functions Vi — Pand vi — ¢ are important in 
correlation work,’ particularly when partial correlations have to be computed from correla- 
tions of the next lower order from formulae such as 


— 
V(1 — — 
The function Vg is useful in work with the binomial distribution, in particular for the 
calculation of standard deviation Vpq/ VN of a percentage p based on a sample of N ob- 
servations. When VNpbq is required, it is often more convenient to multiply g by the integer 
Np (the observed number having a certain characteristic) and use a table of square roots. 
Table II is a small 2-page table of the incomplete -function, known to statisticians 
as the probability integral P of the statistic x (i) (k), and given by 
2(4n)* 
dx. 


Four decimals are given for P, while the parameter m (called the degrees of freedom) takes 


the values 1(1)10, 12, 15, 19, 24 and 30. The argument x/Vn has the range 0-0(0-1)4-1; 
for m greater than about 12 or 15 the interval of 0-1 is too large. 
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Tables III, IV and V. Since no differences are given in Tables I and II, auxiliary tables 
of Lagrangian coefficients are provided for interpolation on a calculating machine. The 
range of argument is 0 to 1; further details are given below. 


Table called points interval decimals 
III Cubic 4 0-001 10 (exact) 
IV — 6 0-01 10 

Vv ptic 8 0-1 11 (exact) 


Unfortunately these tables are bound in with the main tables, so that it is awkward to use 
them for interpolation in Tables I and II. They would have been more convenient under a 
separate cover. Errors are discussed in MTE 24. The accuracy of the various interpolation 
formulae is fully discussed and demonstrated over the whole range of the table. Inverse 
interpolation becomes necessary, as the author wishes to cover cases when (x) or 2(x) are 
required; the processes described are very cumbersome, because the simple methods of in- 
verse interpolation (e), which make use of printed differences, are not applicable here. 

Table VI gives square roots of all integers from 1 to 1000, to 7 decimals up to 99, and 
thereafter to 6; in other words, there are always 8 figures. The layout, in which consecutive 
values are in lines rather than in columns, is not convenient. 

In conclusion it may be said that the value of these tables lies in the accurate record of 
normal deviates at a very fine interval. This should be useful in fundamental work on 
statistical distributions derived from the normal law. The working statistician, on the 
other hand, wants a larger number of statistical functions to fewer decimals, rather than a 
few functions to many decimals. 


Some Related Tables 


(a) W. F. SHEPPARD, The eee Integral (B.A.A.S., Mathematical Tables, v. 7), Cam- 
bridge, University Press, 19 

(b) R. A. FisHeErR and F. Ran ” Shstietiont Tables, second ed., London, Oliver and Boyd, 
1943, Table I. 

(c) FISHER and YartEs, idem, Table II. 

(d) FisHeR and YATEs, idem, Table IX. 

(e) L. J. Comrie, Interpolation and Allied Tables. Reprinted from the Nautical Almanac 
for 1937. London, H. M. Stationery Office, 1936, p. 934-940. 

(f) J. R. Miner, Tables of V1 — r* and 1 — r*. Baltimore, Johns Hopkins, 1922. 

(g) K. PEARSON, Tables for Statisticians and Biometricians, part t th third ed., London, 
Biometrika, 1930, Table I. 

(h) K. PEARSON, idem, Table IT. 

(i) K. PEARSON, idem, Table XII. 

(j) for Statisticians and Biometricians, part 11. London, Biometrika, 
1931, Table 

(k) K. PEARSON, Tables of the Incomplete T-function. London, Biometrika, 1922. 

()) MATHEMATICAL TaBLEs Project, New York, Tables of the Probability Integral, v. 2, 
194 


23 Bedford Square, 
London, W. C. i, England 
H. O. HARTLEY 


1 The quantity x + 5 is equivalent to a table of probits as used in biometric work. 


2 Miner’s table (f) of V¥1 — 1? gives 6 decimals for r = 0-0000(0-0001)1-0000. Its 
layout is less convenient than that of the present table. 


131[K, L]—Epwin Raymonp SMITH (1879- ), “Zeros of the Hermitian 
[sic] polynomials,” Amer. Math. Mo., v. 43, 1936, p. 354-358. 


In Hermite’s original papers, “Sur un nouveau développement en série des fonctions,” 
Institut de France, Acad. Sci., Comptes Rendus, v. 58, 1864, p. 93-100, 266-273, his poly- 
nomials were defined by 

Un(x) = e* 
dx" 


which satisfies the equation U,/" — 2xU,,’ + 2nU, = 0. 


= 
| 
| 
£3 
‘ 
is 
| 
| 
| 
| 
| 
| 
| 


RECENT MATHEMATICAL TABLES 153 


Mr. Smith’s Hermitean polynomials are defined, like those in Jahnke and Emde, 
compare RMT 113, and in W. Hahn's “Bericht’”* by 


H,(x) = (- 


and are of the form 

2 2-4 
= (— v2)"H,.(v2x). 

Mr. Smith states that he had not been able to find in print any list of numerical values of 
the roots of Hermitean polynomials, and that apparently none had been computed except 
for special cases. Hence he computed the zeros for m = [3(1)27; 6D]. 


The functions of the parabolic cylinder, u,(x) are connected with the Hermitean poly- 
nomials by the equation 


un(x) = 
which satisfies the differential equation 
un" (x) + (n + 4 — 32*)u,(x) = 0. 

Since H,(x) and u,(x) vanish for the same values of x, methods due to Sturm are applied 
to this equation for obtaining first approximations to the zeros for small values of ». He 
remarks that while this procedure is perfectly general it is easier for relatively large values 
of m to obtain the approximations by extrapolation from the zeros of polynomials cor- 
responding to m — 1, m — 2, ---, by well-known difference methods. 

He then notes that better values may be obtained from these first estimates by Newton’s 
method of successive approximations. For comparatively small values of m, however, a 
direct application of this method is not feasible; in this case he indicates an adaptation of 
the method suitable for machine calculation. 

CA, 


1W. Hahn, “Bericht iiber die Nullstellen der Laguerreschen und der Hermiteschen 
Polynome,” Deutsche Mathem.-Ver., Jahresb., v. 44, 1934, p. 215-236. It is here pointed 
out, following Hermite, that the H n(x) first occurred in 1836 in a paper of Sturm, and that 
in 1839 they are incidentally mentioned in a paper by Chebyshev. 

J. B. Russell (Jn. Math. Phys., M.1.T., v. 12, 1933, p. 291-297) has defined a Hermite 
function by the relation ¢,(x) = 2%e- 4777 (x) and given a table of ¢.(x) for m = 0(1)11, 
and x = [0.00(0.04) 1.00(0.10)4.00(0.20)7.00, 7.50, 8.00; to at least 5S]. 


132[K, L].—E. R. Smita and ArcuiE HicGpon, “Zeros of the Legendre poly- 
nomials,”’ Iowa State College, Jn. Sci., v. 12, 1938, p. 263-274. Compare 
RMT 92, and 118. 


In this table are given the non-negative zeros of P,,(x), to 6D, from n = 2 to nm = 40 
inclusive. The zeros for n = 2 to m = 7 were abridged from the zeros to 15D in J. W. L. 
Glaisher’s article, B.A.A.S., Report, 1879, p. 53. These zeros may have been found seven 
years earlier. In any case they antedate the values of these same zeros to 16 D, given in 
Hee, Kugelfunktionen, v. 2, Berlin, 1881, p. 15-16, and are in essence equivalent to those 
given by Gauss in 1816. 

In 1934 Mr. HiGpoNn computed the zeros of the polynomials from m = 8 to n = 25 of 
the present article. He used a method based on some theorems by Sturm, referred to in 
RMT 131. Mr. Smith’s computations were for zeros of polynomials from » = 26 tom = 40. 
For this purpose he developed two formulae (p. 263-271) which may be used to compute 
approximate values for the zeros. One of these formulae is intended primarily for the values 
near zero and the other for values near unity, either of the formulae being suitable for 
computing the middle values. The method by which the formulae were obtained is an 
adaptation of the one which was used in a paper by F. ZERNIKE (‘‘Eine asymptotische 
Entwickelung fiir die grésste Nullstelle der Hermiteschen Polynome,” Akad. v. Wetens., 
Proc., v. 34, 1931, p. 673-680). 
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Unit errors in Mr. Higdon’s zero values, from n = 8 to nm = 16, found upon comparing 
them with those of A. N. Lowan, e¢ al., RMT 92, are as follows: n = 8, for .525533, read 
525532; m = 11, for .519095, read .519096; m = 12, for .367832, read .367831; n = 14, for 
-827202, read .827201; » = 16, for .095012, read .095013, and for .281605, read .281604. 

R.C. A. 


133[L].—A. N. Lowan and MILTON ABRAMOWITZ, ‘‘Table[s] of the integrals 
Se Jo(t)dt and fo* Yo(t)dt,” Jn. Math. Phys., M.1.T., v. 22, 1943, p. 2-12. 
17 X 25.2 cm. 


Accompanying these tables is not a single reference to the literature of the subject. 
Integrals of Bessel functions have occurred for a long time in physical investigations and 
until recently adequate tables have been unavailable. In Koppr’s work on the propagation 
of waves along a weighted string! and in Orr’s work on radiation from an alternating 
current,? integrals of Bessel functions of integral order occurred. Estimates of the value for 
large values of the order and argument were given by PoLLaczEK.? Integrals of the Bessel 
functions of the second kind occurred in the work of HAVELocK‘ on ship waves and wave 
resistance. The need of tables was met to some extent by the 7-place tables of Watson® 
which covered the range x = .02(.02)1 and included values at the first 16 zeros of the 
integrand Jo(t) or Yo(#). Four-place tables for the same range were given by GLAZENAP.® 
Graphs of the first integral, divided by x, were given in the papers of DEBYE’ and FiscHER.*® 
This first integral has occurred recently in some work on the intensities of electron dif- 
fraction rings by BLACKMAN® who gives an equivalent definite integral 


RA) = f° sin® [A (1 + + w*) = Jo(2yat 


and a plot of R(A) which indicates that a maximum value 7.2 occurs when A is about 1.25. 
This function R(A) is a particular case of a function of 4 variables studied by BucHHOLz!® 
in his work on the spreading of alternating currents in the earth. 

The tables under review are given for the range x = [0(.01)10; 10D], and are a by- 
product of the computation of Bessel functions of complex arguments. Power series are 
used for the first integral while the computation of the second integral is effected with the 
aid of the first integral multiplied by a logarithm and two power series the first of which 
is evaluated by iterated Taylor expansions, by a method similar to that described in the 
Introduction of Tables of Sine, Cosine and Exponential Integrals, v. 1, prepared by the 
MATHEMATICAL TABLES PROJECT. 


H. B. 


1M. Koppe, Die Ausbreitung einer Erschiitterung an der Wellenmaschine durch einen 
neuen Grenzfall der Besselschen unctionen, Progr. Andreas Realgymn., Berlin, 1899; see 
also Forts. or. v. 30, p. 420-421. 
2W. Mc. F . ORR, “Note on the peeneiten from an alternating circular electric current,” 
Phil. Mag., s. 6, v. 7, 1904, p . 336-34 
3F, PoLLAczeK, “Uber mechanischer Vorgange in einem linearen 
Gitter,” Annalen d. Physik, s. 5, v. 2, 1929, p. 991-1011. 
*T. @. ag ge ‘Studies in wave resistance : the effect of parallel middle body,” 
R. So. London, Proc., v. 108 A, 1925, p. 77-92; ‘Ship waves: their variation with certain 
systematic changes of comm” R. So. ‘London, Proc., v. 136 A, 1932, p. 465-471. 
5G. N. Watson, Theory of Bessel Functions, Cambridge, Univ. Press, 1922, p. 752. 


6S. a GLazeENAP, Matematicheskie i Astronomicheskie Tablitsy, Leningrad, Acad. Sci., 


1932, 
iB Dy Desve, “‘Zerstreuung von Réntgenstrahlen,’’ Annalen d. Physik, s. 4, v. 46, 1915, 

p. 809-823 

8F. A. Fiscuer, “Uber die akustische Strahlungsleistung von ins- 
besondere der Kreis- und Kugelgruppen,” Elek. Nach. Tech., v. 9, 1932, rp. 147-15 

BLackMAN, “On the intensities of electron diffraction rings,” R. So. 
Proc., v. 173A, 1939, P. 68-82. 

10H, BucHHoLz, ‘ ‘Die Wechselstromausbreitung im Erdreich unterhalb einer einseitig 
offenen und unendlich langen vertikalen Leiterschleife im Luftraum,” Archiv f. Elektro- 
technik, v. 30, 1936, p. 1-33. 
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134{L]—ArNoLD Lowan (1898- ), GertTrupDE BLANcH, and 
MiLton ABRAMOWITz (1915— ), “Tables of Jig(x) = /2°J(t)dt/t and 
related functions,” Jn. Math. Phys., M.1.T., v. 22, 1943, p. 51-57. 
17 X 25.2 cm. 


Since this paper contains practically no reference to the literature of the subject, it 
may be well to give indications in this regard, so as to articulate what is new in the tables 
under review. Integrals of the type /-J(x)x*dx were studied in a general way by Lommel! 
but the first use of the integral tabulated in this paper may have been in 1889 when the 
brothers LopGE* obtained the integral 


in some physical investigations. The integral Jio(x) (see Edinb. Math. So., Proc., s. 2, v. 3, 
p. 276) was transformed into a definite integral by LERCH? in 1896. The result involves the 
gamma function of a complex variable and a parameter w whose real part lies between 0 
and 4. The notation adopted here is essentially that of B. van DER Pot‘ who showed that 
the equation 


aJn(x)/an = — n>1 


could be used to calculate the derivative of the Bessel function with respect to n. The 
suffix » was added by P. HumBert® when the Bessel-integral function was defined for an 
arbitrary value of 2 by the equation 


Fine) =— In(Oat/t. 
Humbert gave a generating function for the functions of integral order 
= tFig(x) 


and obtained a recurrence relation 
(nm — 1) Jin-a(x) — + 1)Jings(x) = (2n/x)T a(x). 


He also gave many series and definite integrals for the Bessel-integral functions. In 1938 
Hoac® obtained an expansion for /2°Jo(t)dt/f**! which is derived again by the present 
authors for the special case r = 0. 

The asymptotic expansions which are given in this paper, and in the adjoining one of 
V. G. Situ,’ appear to be new. The coefficients are calculated here to a large number of 
significant figures so that the series can be used for accurate computations. Smith’s ex- 
pansion is convenient for 0 < x < 25 while that of this paper is thought to be better for 
x > 25. 

In Table I the function F(x) = Jio(x) + In(4x) is tabulated with &F(x), 5*F(x) and 
8*F(x), for x = [0(.1)3; 10D]. The function Ji(x) is tabulated with central differences & 
and 8‘ in Table II, for x = [3(.1)10; 10D], and without differences for x = 10(1)22, the 
values for x = 0(.1)3 being given in Table I. Table III contains reduced derivatives of F(x); 
the function A,(x) which is tabulated for m = [0(1)13; 12D], x = 10(1)21, is defined by 
the equations 

A,(x) = F(x)/n!, Ao(x) = F(x). 


These new tables will be much welcomed. 
H. B. 


1E. Lommet, “Uber eine mit den Bessel’schen Funktionen verwandte Funktion,” 
Math. Ann., v. 9, 1876, p. 425-444. 

20. Jj. Lopce and A. Lonce, “The rotation of the plane of paate of light by the 
discharge of a Leyden jar,” Phil. Mag., s. 5, v. 27, 1889, p. 339-346. 
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3M. Lercn [Speculations on some questions in integral equations] (Bohemian), 
Conk Akad., Rozpravy, Class II, v. 5, 1896, No. 23, 16 p.; Forts. Math., v. 27, 1896, p. 233- 

4B. vAN DER Pot, “On the operational solution of linear differential equations and an 
investigation of the properties of these solutions,"’ Phil. Mag., s. 7, v. 8, 1929, p. 861-898. 

5 P. HumBert, “Sur les fonctions de Bessel-intégrales,”’ Inst. d. France Acad. d. Sci., 
Comptes Rendus, v. 195, 1932, p. 854-855; “‘Bessel-integral functions,’’ Edin. Math. So., 
Proc., s. 2, v. 3, 1933, P. 276-285. 

H. A. Hoag, ‘ ‘Equilibrium | of a thin symmetrically loaded, on an 
elastic foundation of infinite depth,” Phil. Mag., s. 7, v. 25, 1938, +7 576-582 

7V. G. Smita, ‘‘An asymptotic expansion of Fete) = Se Jo(t)dt/t,’’ Jn. Math. Phys., 
M.LT., v. 22, 1943, p. 58-59. 


135[L].—A. N. Lowan, and WILLIAM HOoRENSTEIN, “On the function 
H(m, a, x) = exp (— ix)F(m + 1 — ia, 2m + 2; ix).” Jn. Math. Physics, 
M.LT., v. 21, 1942, p. 264-283. 


This paper contains tables of confluent hypergeometric functions for certain complex 
arguments. When normalized in a certain way, the confluent hypergeometric functions are 
also called Whittaker functions; see WHITTAKER and Watson, Modern Analysis, fourth ed., 
Cambridge, 1927, and Amer. ed., New York, 1943, chapter 16. 

In a systematic study of second order ordinary differential equations, the hyper- 
geometric functions play an essential réle. If the differential equation possesses less than 
three regular singular points, the solution can be obtained in elementary functions; with 
three regular singularities the solution, while no longer expressible in terms of elementary 
functions, can be solved by hypergeometric functions. When two of the regular singularities 
coalesce, an irregular singularity is obtained, and the resulting equation has confluent 
hypergeometric functions as solutions. 

In applications such differential equations have arisen in a number of connections. 
H. A. WEBB and J. R. Arrey, in Phil. Mag., s. 6, v. 36, 1918, p. 129-141, suggested a number 
of applications of the confluent hypergeometric function in the theory of elasticity, and in 
problems of flexure and vibration, and published tables and graphs of the function for certain 
real arguments. See MTAC, ‘Mathematical tables in Phil. Mag.,” no. 24. For further 
tables published by Arrry see “‘Mathematical Tables in Reports of the B.A.A.S.,”’ nos. 41, 
43, MTAC, p. 73. Later, R. GRAN OLsson discussed the confluent hypergeometric function 
in connection with the buckling of elastic plates of variable thickness, and the whirling of 
shafts; and included further tables. (Zeit. angew. Math. Mech., v. 16, 1936, p. 347-348; 
Ingenieur-Archiv, v. 8, 1937, p. 81-98, 99-103, 270-275, 373-380.) A collection of curves 
is found in JAHNKE and EmpE, Tables of Functions, New York, 1943, p. 275-282. 

All these applications and tables involved only real arguments and real values of the 
parameters. In the quantum theory of atomic collisions, i.e., in the solution of the 
Schrédinger wave equation for a Coulomb field, the confluent hypergeometric function 
occurs with imaginary arguments and certain complex values of the parameters. (See 
A. SOMMERFELD, Ann. Physik, s. 5, v. 11, 1931, p. 257-330, especially p. 272-273; N. T. 
Morr and H. S. W. Massey, Theory of Atomic Collisions, Oxford, 1933, p. 36.) The present 
tables have been prepared with special reference to these quantum mechanical applications. 

The confluent hypergeometric function y may be defined by the series 


4 alfa +1)---(a+n) xu! 
v(y + 1) 2! + 1)---(y +n) (n+ 1)! 
The corresponding differential equation is 


Fla, y3u) =1+" 


+ (y — u)y’ — ay = 0. 


The tables give seven-place values of H(m, a, x), as defined in the title of the paper, for 
x = 0(1)10 and for a = 0(1)10, m = 0(1)3. 
P. W. KetcHuM 
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136[L].—A. N. Lowan, and H. E. Sauzer, “Table of coefficients in numerical 
integration formulae,” Jn. Math. Phys., M.1.T., v. 22, 1943, p. 49-50. 
17 X 25.2 cm. 


The Bernoulli polynomial B,™(x) of order m and degree m may be defined by means 
of the generating function 


— 1) ™ = /n! 


Its properties are given very fully by NOrLUND,! MiLNE-THomson,? and Davis.? When 
m = n there is an integral representation 


When x = 0 the value of the polynomial is denoted by the symbol B,™ and the numbers 
B, are connected by the recurrence relation 


G)B./s! = (— 1)*B,™ /n! 

Values of these numbers for m = 1(1)5 were given by NORLUND on p. 192 of his paper, and 
on p. 459 of his book there are tables for m = 1(1)12, values being given also for B,‘™/n!. 

On p. 209-210 of his book Davis gives expressions for B,™(0) fer » = 0(1)12 and 
the values of B,‘”) may be deduced by putting m = n. In the present tables both B, (0) 
and B, (1) are given for m = 1(1)20, the values for the higher suffixes being given to twenty 
or more significant figures. The numbers occur as coefficients in Laplace’s formulae of 
numerical integration. The latter set occurs when forward differences are used as in the 
Gregory formula of interpolation, and the former set when backward differences are used. 
The new values provide for much greater accuracy in numerical integration. 
H. B. 


1N. E. NOriunp, “Mémoire sur les polynomes de Bernoulli,’’ Acta Math., v. 43, 1922, 
p. 121-196; Vorlesungen tiber Differenzenrechnung, Berlin, Springer, 1924. 
The Calculus of Finite Differences, London, Macmillan, 1933, 
p- 191-1 
3H. T. Davis, Tables of the Higher Mathematical Functions, Bloomington, Ind., The 
Principia Press, v. 2, 1935, p. 208-210. 


137[L].—J. A. Stratton, P. M. Morse, L. J. Cau and R. A. HUTNER, 
Elliptic Cylinder and Spheroidal Wave Functions, New York, Wiley; 
London, Chapman & Hall, 1941, xii, 127 p. 21 X 27.2 cm. [P. 1-51 are 
an exact reprint of an article by Chu and Stratton in Jn. Math. Physics, 
M.L.T., v. 20, Aug. 1941, p. 259-309. To this were added a title-page 
and foreword (p. i-xii) and “Errata’’ sheet by Morse, and tables, p. 
53-125, by Stratton, Morse, Chu, Hutner, J. B. Fisk, and J. R. Pelham; 
the introduction to the tables, by Hutner, occupies p. 53-76. All the 
extra pages are almost wholly reproductions of typed material.] $1.00. 


The linear differential equation of the second order which, in the normal form, has no 
term involving the first derivative, unity as the coefficient of the second derivative and a 
periodic function P(t) as the coefficient of the dependent variable y, is usually called Mathieu's 
equation, or the equation of the elliptic cylinder, when P(¢) is of the form a + 16g cos 2¢. 
The special case in which a = 1 was used by Laplace in 1777 to iliustrate his method of 
solving differential equations approximately.! He first expanded the solution in powers of g, 
neglecting terms of order higher than g*, and obtained two linearly independent solutions 
involving trigonometrical functions of ¢ some of which were multiplied by ¢ and #. Not 
content with this form of solution he adopted another method ! not unlike the method of 
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the variation of arbitrary constants, and obtained a solution in which trigonometrical func- 
tions of odd multiples of ¢ up to the fifth were multiplied by exponential factors e**¢'. Laplace 
thus introduced the idea of a characteristic index which for his case and for the order of 
approximation adopted was +4g. At this time he was interested in the question of the 
stability of the solar system and his thoughts on this question may have led to the search 
for this second form of his solution. Laplace pointed out, moreover, that his method was 
applicable also to the more general equations with periodic coefficients which occur in 
astronomical investigations. The equation in which P(t) has the form of a cosine series of 
even multiples of ¢ was studied in detail a hundred years later by G. W. Hill and is commonly 
called Hill’s equation because he introduced a novel method of treatment based on the use 
of infinite determinants, a method which was made rigorous by Henri Poincaré and 
von Koch.? 

These linear differential equations with periodic coefficients became important also in 
acoustics and were used by Rayleigh in his paper on the maintenance of vibrations by forces 
of double frequency. Many other writers like A. Stephenson and C. V. Raman have used 
these equations in the study of the stability of a dynamical system or of some peculiarities 
of the oscillations. 

Another very extensive field in which equations of this type occur was opened up when 
in 1837 G. Lamé ‘ focussed attention on the problem of the separation of variables in the 
study of the equation of heat conduction, and in the case of thermal equilibrium studied 
the simple solutions of the partial differential equation having the form of a product of three 
functions, each of a single variable, now called a Lamé product. In 1868 this idea was used 
by E. Mathieu § in the study of the vibrations of an elliptic membrane and it became im- 
portant to find when the differential equation had a periodic solution. Such solutions are now 
known as Mathieu functions and the associated harmonic functions which occur in the solu- 
tion of potential problems are called the functions of the elliptic cylinder. In 1873 Mathieu 
used elliptic coordinates also § for the thermal problem in which there is conduction of heat 
and showed that the separation of variables leads to a more general equation which includes 
the former equation as a special case. This equation is the one on which the theory of sphe- 
roidal wave functions is based. It was studied in some detail by C. Niven * in 1880 while 
Heine said 7 in 1881 that he had been unable to complete the attempt to use in the general 
case the method which proved successful in the treatment of the equation of the elliptic 
cylinder in his first volume. 

Heine succeeded, indeed, in classifying the periodic solutions of Mathieu's equation into 
four classes, the trigonometrical functions in the expansion used being either all sines or all 
cosines and the multiples of ¢ either all even or all odd. This classification has been adopted 
by Whittaker and others as it suggests a convenient notation which, however, is not used by 
the present writers. 

In his early work on Lamé functions Heine used expansions in series of associated 
Legendre functions and in his later work on the functions of the elliptic cylinder he used also 
expansions in series of Bessel functions. This plan was adopted by Niven for the study of 
the spheroidal wave functions and toa slight extent by M. Brillouin * in 1904. Its advantages 
were shown by E. G. C. Poole in 1923 and were particularly emphasized by E. A. Hylleraas ® 
in his lectures in 1937 at the Institute of Henri Poincaré. Not only is the approximation 
more rapid than when power series are used (as in the work of Maclaurin and the associated 
work of Brillouin) but the recurrence relations by which the coefficients are calculated are 
essentially the same. This fact seems to depend on the circumstance that the desired func- 
tions are solutions of a homogeneous integral equation discovered by M. Abraham in 1899 
and in a more general form by E. T. Whittaker,!° who also obtained a corresponding result 
for Lamé functions. 

The associated Legendre functions and Bessel functions used in the expansions are 
limiting forms of the spheroidal wave functions and the use of these expansions is analogous 
to the method of W. Ritz for the direct solution of a problem in the Calculus of Variations. 
It is perhaps on this account, as Hylleraas points out, that the method of computation is 
so effective. 
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This method is essentially that used in the work under review except that the related 
functions of Gegenbauer are used in place of the associated Legendre functions, the notation 
being different from that used by Nielsen.“ A new notation applicable also to the elliptic 
cylinder functions is adopted for the spheroidal wave functions and special values of the 
functions are listed. It is necessary to consider both prolate spheroidal functions and oblate 
spheroidal functions. The functions are also classified according to the type of Gegenbauer 
function or Bessel function used in the expansion. Eight U-functions U,(a, c; z) are defined 
by means of series of different types of Bessel functions and eight V-functions V,(a, c; z) 
are defined by means of different types of series of Gegenbauer functions. Questions of 
convergence are discussed, expansions are obtained in the form of power series and relations 
are obtained among the U-functions, among the V-functions and between the U-functions 
and V-functions. The corresponding solutions of the wave-equation are denoted by an 
S-notation and an R-notation. Different kinds of functions are distinguished which in turn 
are classified as even or odd. Combinations are formed also to correspond to the Hankel 
functions and serve for the solution of radiation problems in electromagnetic theory. Func- 
tions are also described as angular or radial, the angular functions being periodic in ¢ where 
z = cost, The separation constants are calculated with the aid of continued fractions by a 
method due to Lindstedt, which was used successfully by Hough in 1897 in an application 
of harmonic analysis to the dynamical theory of the tides. This method was recommended 
by Poole in 1923 and was used successfully by Ince and Goldstein in their work on the 
tabulation of Mathieu functions from 1927 on. 

Another change of notation is effected by the introduction of a normalization factor, 
and the coefficients D,!' which are tabulated to 5D in the first 5 tables for / = 0(1)4 are the 
coefficients in the expansion in a series of Bessel functions J,(cz) of a function Jer(c, z) which 
is a multiple of the U-function of the Mathieu type. Seven coefficients are tabulated and 
values are given also for the separation constant 5;. The next 4 tables give for / = 1(1)4 the 
values of the coefficients in the expansion in a series of Bessel functions J,(cz) of a function 
Joi(c, z) which is a multiple of the U3(c, z) function of the Mathieu type. Seven coefficients 
and the separation constants are given. In the list of errata a misprint is noted in the equa- 
tion (310). The values of the separation constants supplement those of Goldstein and Ince 
as they are for small values of ¢ and small interval while c is used as variable instead of a 
multiple of c*. 

The next table gives the separation constant A, as in the equation on p. 62, for the pro- 
late spheroidal functions, to 5D, for c = 0(.1).6(.2)5. The meaning of the suffixes is not 
given; one of them is presumably m, the other /. The values m = 0, 1, 3 are chosen, / being 
zero in the last case and ranging from 0 to 3 in the first case and from 0 to 2 in the second. 
This meaning is more or less confirmed by the next table of the coefficients d in the expan- 
sions of the S- and R-functions for the prolate spheroid. The range of c is the same as before 
but m also has the value 2 which is missing in the table of the A’s. About 16 coefficients are 
given for each case. The remaining tables refer to the oblate spheroidal wave functions. The 
table of B’s on p. 107 refers to the B in the equations on p. 69, and it i- to be noted that 
the value m = 2 is now given in the table but only with / = 1, whereas in the table of the 
coefficients f the value / = 0 is also considered. 

The tables are evidently intended to be used chiefly for the solution of electromagnetic 
problems. In the problems of wave mechanics connected with the hydrogen molecule ion 
the internal equation is the same as the equation considered in the present paper but the 
external equation is a little different. The short table calculated by Hylleraas is consequently 
not comparable with these tables. The only comparison that can be made is with the few 
numerical results obtained by Maclaurin ® and M. Brillouin * in their work on electrical 
oscillations and with the new results obtained recently by Page and Adams." It can, then, 
be said that these tables form a distinct addition to our knowledge of spheroidal wave 
functions. 

H. B. 


1P. S. pE Lapace, “Mémoire sur l’intégration des équations différentielles par approxi- 
i ’ Mém. Acad. roy. sci., Paris, 1780, or Oeuvres, v. 9, p. 357-380, especially p. 359-362, 
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on A :* —_ History of the American Mathematical Society 1888-1938, New Y ork, 
1938, p. 120. 

HELGE von Kocu, “Sur une application des déterminants infinis 4 la théorie des équa- 
tions differentielles linéaires,’’ Acta Math., v. 15, 1891, p. 53-63; ‘‘Sur les déterminants 
infinis et les équations differentielles linéaires,’’ Acta Math., v. 16, 1892-3, p. 217-295. 

3 RayLeiGH, “On maintained vibrations,’’ Phil. Mag., s. 4, v. 15, 1883, p. 229-235, or 
Scientific Papers, v. 2, 1900, p. 188-193; ‘‘On the maintenance of vibrations by forces of 
double frequency, and on the propagation of waves through a medium endowed with a 
oo — Phil. Mag., s. 4, v. 24, 1887, p. 145-159, or Scientific Papers, v. 3, 
1902, p. 1-14. 

‘ & Lamé&, ‘Mémoire sur les surfaces isothermes dans les corps solides homogénes en 
équilibre de température,”’ Jn. de Math., s. 1, v. 2, 1837, p. 147-188; ‘Mémoire sur l’équilibre 
des températures, dans les corps solides homogénes . . .,”’ Jn. de Math., s. 1, v. 4, 1839, 
p. 126-163, 351-385. 

5 E. L. Martuieu, “Mémoire sur le mouvement vibratoire d’une membrane de forme 
mg Jn. de Math., s. 2, v. 13, 1868, p. 137-203; Cours de Physique Mathématique, 

aris, 1873. 

6 C. Niven, “On the conduction of heat in ellipsoids of revolution,’’ R. So. London, 
Trans., v. 171, 1880, p. 117-151. 

7 HEINE, Handbuch der Kugelfunktionen, second ed., Berlin, v. 2, 1881. 
one ” M. BRILLouIn, Propagation de lV’Electricité, Paris, Hermann, 1904, Ch. 6, p. 376, 

, etc. 

*E. A. Hytieraas, “Equation d’ondes d’un électron dans le champ de forces de deux 
noyaux atomiques, probléme de l’ion moléculaire d’hydrogéne,”’ Inst. Henri Poincaré, 
Annales, v. 7, 1937, p. 121-153. 

1 E. T. WHITTAKER, “‘On the functions associated with the elliptic cylinder in harmonic 
analysis,’ Intern. Congress Math., Cambridge, 1912, v. 1, p. 366-371. ‘On Lamé’s differ- 
ential equation and ellipsoidal harmonics,” London Math. So., Proc., s. 2, v. 14, 1915, 
p. 260-268 or Modern Analysis, p. 564; ‘‘On an integral equation whose solutions are the 
functions of Lamé,’’ R. So. Edinburgh, Proc., v. 35, 1915, p. 70-77. See also J. H. PRIESTLEY, 
“On some solutions of the wave equation,” London Math. So., Proc., s. 2, v. 20, 1922, p. 
37-50. J. L. SHarma, Jn. de Math., s. 9, v. 16, 1937, p. 199-203, 355-360. 

11 N, NIELSEN, Théorie des Fonctions Metasphériques, Paris, 1911. 

2 R. C. Macrtaurin, “On the solutions of the equation (y? + k)y = 0 in elliptic co- 
ordinates and their physical applications,” Cambridge Phil. So., Trans., v. 17, 1898, p. 
41-108. For numerical results see p. 92-95, 98, 106, 108. 

13. Pace and N. I. Apams, Jr., “The electrical oscillations of a prolate spheroid, 
part I,”’ Phys. Rev., s. 2, v. 53, 1938, p. 819-831. 
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References have been made to Errata in RMT 121 and 122 (War 
Dept.), 123 (BRowN & SHARPE), 125 (SERVICE GEOGRAPHIQUE), 132 
(Hicpon), and in 118 (HOLLAND, JONEs and LAMB; BYERLY); see also in the 
second article of this issue, nos. 8 and 45. For errors made by MEIssEL and 
WaTtTSON see the first article in this issue. 

21. J. BourGet, {Tables of the first nine roots of J,.(xn) = 0, s = [0(1)5; 3D]}, 

Paris, Ecole Normale Sup., Annales, v. 3, 1866, p. 82-87. 


On comparing the 54 entries of these tables to 3D, with the tables to 5S or 6S for 
s = [0(1)10(5)20(10)50 --- 1000], of J. R. Arrey, B.A.A.S., Report, 1922, p. 271 (see 
MTAC, p. 72), it was found that 27 of the Bourget entries were erroneous. Four of the 
worst errors were as follows: J3(x»), for 32.050, read 32.065; J4(xs), for 33.512, read 33.537; 
Js(x1), for 8.780, read 8.772; Js(xs), for 34.983 read 34.989. All 27 of these errors are faith- 
fully reproduced in each of the five editions of JAHNKE and Empe, Table of Functions; see 
RMT 113. So also for Table V (p. 302) in Gray and MATHEWs, Treatise on Bessel Functions, 
second ed., London, 1922; for Table XXXIII (p. 82) in J. B. DALE, Five-Figure Tables of 
Mathematical Functions, London, 1903; and for Table V (p. 286) in W. E. Byerty, An 
Elementary Treatise on Fourier’s Series . . . , Boston, 1895. These 27 errors (and one 
more added) are reproduced with equal faithfulness in RAYLEIGH, The Theory of Sound, 
v. 1, London, 1877, p. 274 (also in the German ed., v. 1, Brunswick, 1879, p. 364, and in 


= 
| 
: 
a | 
) 


MATHEMATICAL TABLES—ERRATA 161 


various other English prints such as v. 1, 1929, p. 330.) Compare J. R. Arrey, Phil. Mag., 
s. 6, v. 32, 1916, p. 10. F. Bowman, Bessel Functions (1938) also copied Bourget errors. 
GA. 


22. [L. J. ComriE], Seven-Figure Trigonometrical Tables for every Second of 
Time, London, 1939. Compare MTAC, p. 43. 


In addition to the errors already included in the Errata Slip are the following: P. 74, 
d for tan 2" 05", for 995-1000, read 995-1001 ; d for cot 3" 54™, for 1000-995, read 1001-995. 
SUPERINTENDENT of H. M. Nautical Almanac Office 


23. P. R. E. JAHNKE and F. Empe, Tables of Functions, first ed., 1909. 
Compare RMT 113. 


P. 54, F(35°, 30°), for 0.6220, read 0.6200. P. 137, table expressing Jo(xVi) in the form 
a + ib, there are 10 errors of sign in the imaginary parts from x = 5.1 to x = 6.0 inclusive 
(F. F. P. Brsacre, Phil. Mag., s. 6, 1923, v. 45, p. 1035.) 


24. T. L. KeLtey, The Kelley Statistical Tables, 1938. Compare RMT 130. 


The Tables III, IV and V (of 4-point, 6-point and 8-point Lagrangian coefficients) have 
been compared with tables of these coefficients that I made about 12 years ago on a Bur- 
roughs machine. In the 2000 values of Table III there is only one error: 

Page 121, p = 0.155, col. Co, for .90033 79375, read .90033 69375. As the values are 
exact, the question of rounding-off errors does not arise. In Table IV there are no errors, 
when we take into account the fact that in 15 cases the rounding off has been (quite properly) 
changed by a unit in order that the sum of the coefficients shall be 1. The 72 values in 
Table V are all exact and correct. 


L. 


25. MATHEMATICAL TABLES Project, New York, I. Tables of the Exponential 
Function e*, 1939; I1. Tables of Circular and Hyperbolic Sines and Cosines 
for Radian Arguments, 1940; compare RMT 89; III. Tables of Sine, 
Cosine and Exponential Integrals, v. 1, 1940; IV. Table of Natural Loga- 
rithms, v. 4, 1941; V. Tables of Probability Functions, v. 1, 1941; compare 
RMT 91 and MTE 16; VI. Table of Arctan x, 1942; compare RMT 90. 
I. P. 202, heading, for ‘‘descending,’’ read ‘‘ascending’’; p. 282, heading, for o*, read e~*; 
p. 335, argument, read .5486. 
II. P. 358, argument, read 1.7850. 
III. P. xv, footnote, for “result in error by six units,” read “result in error by two units”; 
also, for (1/24)p(p — 1)(p — 2)(b — 3), read (1/24) p(p* — 1)(p — 2). 
IV. P. 440, argument, read 9.3890. 
V. P. 131, argument, read 0.6496. 
VI. P. 36, arguments, read 3.456, and p. 71, read 6.905. 
A. N. Lowan 


26. C. W. MERRIFIELD, “‘The sums of the series of the reciprocals of the 
prime numbers and of their powers,” R. So. London, Proc., v. 33, 1881, 
p. 4-10. 


Mr. K. W. Miller, of Chicago, has called my attention to the following errata (p. 10): 


n Sor read 
3 17096 17476 
5 50164 50174 

11 69104 69105 


13 26973 26983 
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There is also a last figure error in Merrifield’s preceding table (p. 9) of log S, for m = 13 
where 1 should be 2. All these errata are in Gram’s reproduction of these tables [K. Danske 
Videns. Selskabs, Skrifter, s. 6, Naturd. og Math., v. 2, p. 269]. Their discovery is the result 
of collation with Davis’ excellent tables of these sums to 24 decimals in his Tables of the 
Higher Mathematical Functions, Bloomington, Ind., v. 2, 1935, p. 249-250. The verification 
mentioned by Glaisher (/.c.) has been applied to the tables of Davis and Glaisher by the 
writer without uncovering any error. 

D. H.L. 


1 The first error was discovered by J. W. L. Glaisher, Quart. Jn. Math., v. 25, 1891, 
p. 0. and listed in the writer’s Guide to Tables in the Theory of Numbers (without due 
it). 


27. J. T. Peters, Sechsstellige Tafel der trigonometrischen Funktionen,.. . , 
Berlin, 1929. 


In acknowledging the list of errors in this work which I had sent to Peters (see MTE 14), 
he wrote on 23 May 1939 that a new proof had been read for the second edition (which 
appeared later in 1939), and that this had confirmed his suspicion that Herr Schoch, who 
is mentioned in the Introduction as one of the readers, had not been sufficiently painstaking 
with his share (except in the case of the last three figures of each function), and that many 
errors remained. I felt that this warning should be passed on to anyone who may be using 
the first edition. 


28. HERMANN BRANDENBURG, Sechsstellige trigonometrische Tafel, Leipzig, 
Lorentz, 1932. xxiv, 304 p. 19.9 K 28.3 cm. 


This table gives 6-figure natural values of the four principal trigonometric functions at 
interval 10’, with an auxiliary table of cotangents up to 3° at interval 1’. Like the similar 
7-figure tables by the same author (Leipzig, Lorentz, 1923 and 1931), it is based on the 
Opus Palatinum of Rheticus (1596). It was compared soon after publication with Andoyer’s 
Nouvelles Tables Trigonométriques Fondamentales (Paris, Hermann, 1915 and 1916). The 
early cotangents were compared with the values computed by the undersigned for Peters’ 
Achtstellige Tafel der trigonometrischen Funktionen fiir jede Sexagesimalsekunde des Quad- 
ranten (Berlin, Reichsamt fiir Landesaufnahme, 1939; Amer. ed., Ann Arbor, Mich., 1943). 
These comparisons and other examinations led to the discovery of the errors listed below, 
only the first of which is serious. This work is, therefore, of a very high standard of accuracy, 
whereas the first edition of the 7-figure tables (1923) was full of errors. Brandenburg, how- 
ever, by offering rewards for any errors discovered, and in other ways, laboured incessantly 
to remove this blemish. The second edition of the 7-figure table (1931) was almost free from 
error, and this volume has fewer errors than the corresponding volume by Peters of all six 
functions to six figures at the same interval (see MTE 14, 27). 


page angle, etc. for read 
4 cot 0°16'24” 208,6171 209,6171 

95 Argument col. for P.P. 1340 3 (line 2) 2 
150 tan 20°09'10” Delete asterisk 
193 Arguments 27°18'40” and 50” To be interchanged 
225 P.P. for 7/10 of 168 117,7 117,6 
302 Beispiel, line 3 mau man 
300 p° 170 32 170 33 
300 p’ 2192 2199 
300 p” 152 197 


The errors in p°, p’ and p” are due to multiplying 1/x (which is correct) by 180, 10,800 
and 648,000, and writing too many significant figures. There are also end-figure errors in 
(}x)*, especially for the higher values of m. The value of e on the same page is, of course, 
correct to 19 decimals only, since the largest neglected term begins in the 20th decimal. 
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UNPUBLISHED MATHEMATICAL TABLES 


We have also referred to unpublished mathematical tables of (a) Si- 
BERSTEIN, in the second article of this issue; of (b) PeETERs, in RMT 125, 
and in N 13; and of (c) PETERS and ComriE in RMT 127. 


15[A].—SoLomon ACHILLOVICH JOFFE (1868-— ) Table of the first 100 fac- 
torials, and the method of its computation. Ms., prepared in January 1921, 
in possession of the author; in 1921 photographic copies were placed in 
the libraries of the Edinburgh Mathematical Society, the University of 
Edinburgh, and the Royal Society of Edinburgh, and in 1943, in the 
Library of Brown University. Compare UMT 10. 


A table of the first 50 factorials is contained in a paper by FRANK Rossins, “Factorials 
and allied products with their logarithms,” R. So. Edinb., Trans., v. 52, part 1, 1917, p. 
167-174. This factorial table was extended by the writer shortly after it was brought to 
his attention. Beginning with m = 51, the multiplications (and additions for verification) 
were made on the Electric Ensign Machine. Each value m! was obtained from (nm — 1)! 
through multiplication by ”, the multiplicand being broken up into 9-figure groups. The 
resulting product was verified by the congruence ! =0 (mod 10° — 1). E.g. for 68! 
2480 + 035542 + 436830 + 599600 + 990418 + 569171 + 581047 + 399201 + 355367 
+ 672371 + 710738 + 018221 + 445712 + 183296 = 5999994. 

Each quinary value n!, for n = 5k, was computed a second time from (m — 5)! through 
multiplication by 2, the latter values having been obtained by successive addition from 
the table of 0, 5, 10, 15, 20 and 25). In multiplying (n — 5)! by n™, the latter 
was taken as the multiplicand and (m — 5)! was separated into 6-figure groups for multi- 
pliers. The final value 100! was also verified by the congruence 100! =— 1 (mod. 101). 

S. A. JorFEe 

515 West 110th St., 

New York City 


16{A, E, L].— Miscellaneous Tables. Mss. prepared by, and in possession of, 
the MATHEMATICAL TABLES Project, 50 Church St., New York City. 
None of the items listed below is intended for publication. 


1. Square roots of the integers up to 100 to about 18S. 

2. The functions x*/n! and x*/n(m!) for x = [0.00(0.01)2; 15D], and m = 1, 2---p, 
where p depends on x and is such that the functions in question are of the order of mag- 
nitude of 

3. Values of x"/n! for x = [1(1)10; 10S], and m = 1(1)40. 

4. The first fourteen derivatives of the function (1/ V2n)e-¥*, for x = [0(0.1)8.4; 20D]. 


5. The first fourteen derivatives of the function (2/ Vx)e=* for x = [0.0(0.1)8.0; 20D}. 

6. The first ten derivatives of the sine and cosine integrals for x = [10(1)100; 12D]. 

7. The first eighteen derivatives with respect to r of the functions Jo(z) and J;(z), 
where z = re”, for @ = 0°(5°)90° and r = 0.0(0.1)10. The derivatives are given to a number 
of places (varying from 8 to 15) adequate for the subsequent subtabulation of the functions 
to about 12D. 

8. The first eighteen derivatives of the functions x~’J,(x) and x~*J,(x) for 
x = 0.0(0.1)10 and » = + 3, + 3, +}, + }. These derivatives are given to a number of 
places (varying from 8 to 15) adequate for the subsequent subtabulation of the functions 
to about 12D. 

9. The first fifteen derivatives of J,(x) and I,(x) for x = 10.0(0.1)25. These derivatives 
are given to a number of places (varying from 8 to 15) adequate for the subsequent sub- 
tabulation of the functions to about 12D or 12S. 
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10. The first eighteen derivatives with respect to r of the functions 


So(z) = Yo(z) — (2/)Jo(z)(Inz + vy), 

Si(z) = Vi(z) — (2/4) Ji(z)(Inz + y) + 2/x2, 
where z = re”, for @ = 0°(5°)90° and r = 0.0(0.1)10. These derivatives are given to a 
number of places (varying from 8 to 15) adequate for subsequent subtabulation of the 
functions to 12D. 

11. Modulus and argument of Jo(z) and Ji(z), where z = re”, for @ = 0°(5°)90° and 
r = [0.0(0.1)10; 10D or 10S}. 
A. N. Lowan 


17[B].—H. T. Davis, Tables of powers. Mss. in possession of the author. 
Compare Q 5. 


I. Tables of reciprocals, values of 1/n, for m = [100(1)1000; 17D], together with dif- 
ferences. The following work giving results to 7D is well known: W. H. Oakes, Table of 
the Reciprocals of Numbers from 1 to 100,000 . . ., London, 1865; there is also the more 
recent work, Barlow's Tables . . . , fourth ed. by L. J. Comrie, London, 1941 (RMT 82) 
which gives the reciprocals of all integers up to 12500, to 7S, with differences; “‘Reciprocals 
of primes less than 10,000,” in units of the tenth and eleventh decimal, are given with dif- 
ferences in B.A.A.S., Mathematical Tables, v. 5, Factor Table, London, 1935, p. 284-291. 
II. Tables of reciprocals of squares and cubes, values of 1/n? and 1/n*, for m = [100(1)1000; 
16S-18S]. There are tables for » = [50(1)200; 10D-13D] in K. PEARson, Tables for Statisti- 
cians and Biometricians, part 2, London, 1931, p. 244-245. III. Tables of reciprocal fourth and 
fifth powers, values of 1/n‘* and 1/n5, for n = [100(1)1000; 18S]. The values of 1/n?, to 32D, 
for m = 2(1)100, p = 1(1)16, are given in J. T. PETERS, Zehnstellige Logarithmentafel, v. 1, 
Berlin, 1922, Anhang, by Peters and Stein, p. 36-57. IV. Tables of the functions x™(1 — x*)" 
for x = (0.00(0.01)1.00; 14D]; m = 0, 1, and 2; m = 1/2, 3/2, 5/2, 7/2, and 9/2. These 
tables were used in computing the values of the integrals of the Legendre associated func- 
tions. In his Tables of (1 — r*)* and 1 — r* for use in Partial Correlation and Trigonometry 
(Baltimore, Md., 1922), J. R. MINER gives the values of (1 — r*)! for 


r = [0.0001(0.0001).9999; mostly to 6D]. 
H. T. D. and R. C. A. 


18[L].—Tables of the Associated Legendre Functions and their First Derivatives, 
Mss. prepared by, and in possession of, the MATHEMATICAL TABLES 
Project, 50 Church St., New York City. 


These are six-place tables of the functions P,”(x) and Q,(x) and their first derivatives 
for real and purely imaginary arguments, integral and half-integral values of m, and integral 
values of m. Details of the tabulation are as follows: 

Tables of P,"(cos 6) and dP,™(cos 6) /d0 for @ = 0°(1°)90°; m = 0(1)10; m = 0(1)4; m =n. 
Tables of P(x), dPn™(x)/dx, Qn™(x), dQn™(x)/dx, for x = 1.0(0.1)10, = 0(1)10; 

m = 0(1)4;m=n. 

Tables of i-"P,™(ix), i-"dP,™(ix)/dx, Qx™(ix), and dQ,™(ix)/dx, for x = 1.0(0.1)10; 

n = 0(1)10; 10; m = 0(1)4; m =n. 

Tables of and dQn4,(x)/dx, for x = 1.0(0.1)10; m = 0(1)4; 

m = 0(1)4 (even where m is greater than 7). 

Comparison of the contents of these manuscripts with the tables referred to in H. T. 
Davis’s review of Mursi’s table of P,™(x) (RMT 118) reveals the following facts: 
(a) Tallqvist’s tables of P,.™(cos 0) cover a greater range in m and give more decimals than 
ours. (b) The function dP,,”(cos @)/d@ is not given by Tallqvist; (c) Mursi’s tables of P»™(x) 
for x less than unity, are completely complementary to our tables for real x. (d) There is 
some slight overlapping between our values of dP,™(cos 6)/d@ and those in the Tallqvist 
manuscript. 

A. N. Lowan 
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19[L].—Table of the Bessel Functions Ko(x) and K,(x) for Small Arguments. 
Ms. prepared by, and in possession of the MATHEMATICAL TABLES 
Project, 50 Church St., New York City. 


The B. A. A. S. published volume of tables of Bessel Functions contains a table of 
Ko(x) and Ki(x) for x = 0.00(0.01)5. Eight significant figures are given and second dif- 
ferences are tabulated for x > 0.3. For x < 0.3 the auxiliary functions E(x), E:(x), Fo(x) 
and F,(x) are also tabulated for purpose of interpolation; the last mentioned functions 
being defined from the formulae 


Ko(x) = Eo(x) + Fo(x) logx, Ki(x) = E,(x)/x + Fi(x) logx. 


Because of the frequent use of these functions in various branches of research, the Project 
has subtabulated the functions above mentioned in the range from 0 to 1. The ranges and 
intervals for the various functions are as follows: For Ko(x) and Ki(x), x = 0.0000(0.0001)- 
0.0300(0.001)1.000. For Eo(x), E:(x), Fo(x) and Fi(x), x = 0.000(0.001)0.030. All entries 
are given to seven significant figures with first and second advancing differences. 

A. N. Lowan 


20.[L].—Bessel Functions of the Second Kind for Complex Arguments. Mss. 
prepared by, and in possession of, the MATHEMATICAL TABLES PROJECT, 
50 Church St., New York City. Film copy at Brown University. 


These are ten-place tables of Bessel functions of the second kind, Yo(z) and Y;,(z), 
for complex arguments z = re”, where r = 0.00(0.01)10, @ = 0°(5°)90°. 

Bessel functions of orders 0 and 1 of the first and second kinds are encountered in the 
general solution of boundary value problems arising in the theory of potential, heat con- 
duction, and wave motion, when the domain is a cylinder or a cylindrical annulus. They 
occur in particular in the problem of propagation of electromagnetic waves, the theory of 
the skin effect for poorly conducting wires, the oscillatory motion of a sphere in a viscous 
medium, the vibration of a heavy chain in a resisting medium, the determination of the 
frictional and gliding coefficients of a fluid in rotary vibrational motion and in many other 
boundary value problems. They have further theoretical importance in connection with 
the problem of conformal mapping and the evaluation of contour integrals involving Bessel 
functions. 

No tables of Bessel functions of the second kind for complex arguments are to be 
found in the literature except for arguments along the 45° ray, which are known as ker and 
kei functions. These functions (or functions closely related to them) have been tabulated 
by H. G. SavincE, Phil. Mag., s. 6, v. 19, 1910 for x = [1(1)30; 4S]; by F. T6LKe, Besselsche 
und Hankelsche Zylinderfunktionen nullter bis dritter Ordnung vom Argument rvi, Stuttgart, 
1936, for r = [0.00(0.01)21.00; 4S]; by H. B. Dwicut, Tables of Integrals and Other Mathe- 
matical Data, 1934, for x = [0.0(0.1)10; 7D-9D]; by N. W. McLacuian and A. L. MEYERs, 
Phil. Mag., s. 7, v. 18, 1934, in polar form, for z = [0.0(0.1)10; 5D to 8D]; by K. Havasu, 
Fiinfstellige Funktionentafeln, Berlin, 1930, p. 111-118. Shorter tables are given by N. W. 
McLacuian, Bessel Functions for Engineers, Oxford, 1934; by E. JAHNKE and F. EmpE, 
Tables of Functions with Formulae and Curves, New York, 1943; and by B.A.A.S., Report, 
1916. 

A. N. Lowan 


MECHANICAL AIDS TO COMPUTATION 


8[Z].—S. L. Brown and L. L. WHEELER, ‘‘Use of the mechanical multi- 
harmonigraph for graphing types of functions and for solution of pairs 
of non-linear simultaneous equations,’’ Review of Scientific Instruments, 
v. 13, 1942, p. 493-495. 
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This note indicates the application of the Brown harmonic analyzer (described in 
MTAC, p. 128-129) to the drawing of “‘complex Lissajous figures” and their application. 
Such figures are curves whose parametric equations express x and y as trigonometric poly- 
nomials in 6, such as 


x = cosi3@ + cos158, y =— sini3é + sini5@, 


which is clearly a rose of 28 petals. Nine drawings, some of them quite complicated, are 
reproduced in this article. 
If f(x, y) and g(x, y) are polynomials, the machine can be used to solve the system 


f(x, = h, g(x, 9) =k, 
for x and y by setting 


(A) x = y =rsiné@ 


and thereby transforming f and g into a pair of functions of 7 and 0, which, for a fixed (trial) 
value of 7, become trigonometric polynomials. If the corresponding complex Lissajous 
figure passes through the point (h, k) for a certain value of 6, this value, and the value selected 
for r, determine x and y by (A). If not, a new trial value of 7 is interpolated. 

D. L. 


9[Z}.—-S. L. Brown and L. L. WHEELER, ‘‘The use of a mechanical syn- 
thesizer to solve trigonometric and certain types of transcendental 
equations, and for the double summations involved in Patterson con- 
tours,” Jn. Applied Phys., v. 14, 1943, p. 30-36. 


The methods and results of applying the Brown analyzer (described in MTAC, p. 
128-129) to the problem of solving such equations as 


x + sinx = 3 
log (2x? + Sx? + 10x + 15) = 5/x 
3 tan’x + 2 cos** + sinx + cscx = 5 


are discussed in detail in this paper. In each case a trigonometric polynomial is obtained 
by some device and the machine used to graph a curve from which the desired roots are 
read off. In the case of the first equation the left member is expanded in powers of x and 
then x is replaced by 2.5 cos@. Powers of the cosine are then reduced to trigonometric 
polynomials. 

Two beautiful graphs are given of Legendre’s P,(cos0) (m = 1, ---, 11), and certain 
associated functions P,™)(cos@). Unfortunately no scales are indicated. Finally there is a 
discussion of the double Fourier series 


n m 
A(x,y) = 2 an cos2ehx cos2aky 
k=0 h=0 


used in the Patterson method for determining interatomic distances in crystals. Here the 
machine is set to draw a series of m + 1 curves 


= LDancos2xhx. (k = 0, ---,n). 


h=0 


Then, reading 5 values of C,z on each curve (for x = 0(.02)1), the machine is set to draw 
the 5 curves 


A(x, y) = Cyz cos2aky. 


These curves can then be used to prepare a contour map of the function A(x, y). 
D. 8.1. 


| 
| 
| 
= 
! 
| 


ly- 


1) 


l= 


(7 


NOTES 167 


10[Z].—-R. L. Dretzoxp, “‘The isograph a mechanical root fincer,”” and R. O. 
MERCNER, “The mechanism of the isograph,”” Bell Laboratories Record, 
v. 16, 1937, p. 130-134 and p. 135-140. 


These short articles describe the theory, practice and construction of a large instrument 
designed in 1928 by T. C. Fry primarily for the solution of algebraic equations of degrees 
= 10. The theory behind the machine is briefly set forth in a review of the S. L. Brown 
harmonic analyzer, to which review the reader is referred (MTAC, p. 128-129). The 
mechanical methods of generating and adding the ten pairs of motions a,r* cos 76, 
a,r" sin n@ and the drawing of this sum are practically the same as those of the Brown 
analyzer. There is one difference however. In the latter machine the two motions of any 
one frequency may be set with any preassigned phase difference, while in the Isograph 
this difference is inexorably 90°. This means that the use of the Isograph is restricted to 
equations with real coefficients. An interesting feature of the Isograph is a set of ten 
cylindrical slide rules driven with frequencies from 1 to 10, for computing the values of 
anr” (n = 1, ---, 10) which are the ten amplitudes to be set on the cross heads. These 
amplitudes must not exceed an inch and a half (as compared with over 6 inches in the Brown 
analyzer) so that frequent changes in scale must be required for accommodating rather 
different values of r. 

The photographs and diagrams illustrating these articles have been put on 12 slides, 
and a film has been made. Since these may be borrowed from the Bell Laboratories (New 
York), together with mimeographed material about the machine, there is the basis for an 
interesting lecture about the Isograph for mathematics classes and club meetings. 

D. H.L. 


NOTES 


11. ALFRED LopGe.—In referring to sketches of two tablemakers in 
Who was Who 1929-1940, London, 1941 (N 1) a third one was overlooked, 
namely: of ALFRED LopGE (1854-1937), brother of Sir Oliver Lodge (1851-— 
1940), long principal of the University of Birmingham, and of Sir Richard 
Lodge (1855-1936), professor of history at the University of Edinburgh. 
Alfred Lodge was Fereday fellow of St. John’s College, Oxford, 1876-91, 
professor of pure mathematics at the Royal Indian Engineering College, 
Cooper’s Hill, 1884-1904, and assistant master (mathematics) of Charter- 
house from 1904 until his retirement in 1919, at the age of 65. He was 
president of the Mathematical Association, 1897-98, and scores of his 
communications appeared in the Mathematical Gazette. For some years he 
was secretary and recorder of Section A of the B.A.A.S. He was a member 
of its Committee on Calculation of Mathematical Tables for nearly fifty 
years; except for Airey, more of his tables were published than those of any 
other member (see MTAC, p. 75). He was a joint author of the Committee’s 
Factor Table (1935), of which three independent copies were made; compare 
N 13. Since Lodge alone had been on the Committee from the day on which 
Bessel functions were first discussed up to 1937, when the Committee pub- 
lished its first volume of tables of Bessel Functions, the volume was dedicated 
to him. It appeared, however, a few days after his death. See also The Times, 
London, 6 Dec. 1937, p. 14; and C. O. Tuckey, “Alfred Lodge,”’ Math. 
Gazette, v. 22, 1938, p. 3-4. 


12. INTERPOLATION.—If there is suspicion of the accuracy of an entry 
in a table with equal argument-intervals, the suspicious entry %#) may be 
checked in terms of the six adjacent entries u_3, u—2, u_; and 1, 42, %3, by the 


4 
in 
| 4 
re 
Is 
d 
il 
) 
1 
L 


168 NOTES 


formula uo = 0.75(u_1 + u1) — 0.30(u_2 + u2) + 0.05(u_s + us) — A®u_;/20. 
If sixth differences are negligible, as is often the case in published tables, 
the term A‘*u_3/20 is omitted, as a zero remainder-term. 

In many tables an earlier difference than the sixth may be neglected, yet 
the formula above is to be preferred to a simpler one using less than six 
entries, because no division on the machine is needed. [This is because the 
coefficient of the central term in (E — 1)® is 20; for no other power has this 
coefficient powers of 2 and 5 as its only factors]. 

This formula is fairly obvious when pointed out, but I have never seen 
it in print. 

C. R. CosENs 
Engineering Laboratory, 
Cambridge Univ., England 


13. JOHANN THEODOR PETERS.—An obituary notice by W. R. A. Klose 
in Z. angew. Math. u. Mech., v. 22, Apr. 1942, p. 220, states that ‘Jean 
Peters,” the greatest mathematical table-maker of all time, died on 24 
August 1941. Since he was born on 25 August 1869 his span of life com- 
pleted almost exactly 72 years. His full name, as given above, is as he sent 
it for publication in Poggendorff’s Biographisch-literarisches Handwérterbuch, 
v. 6, part 3, 1938. In his published mathematical tables he always gave his 
name as “Dr. J. Peters” or ‘‘Prof. Dr. J. Peters.’’ In the Portrdtgallerie, 1931 
(MTAC, p. 12) the name is ‘Johannes Peters.”” Since we know the name 
which Peters himself regarded as complete for the last decade at least, 
whence the name ‘“‘Jean’’? 

Peters was born in the rather genial atmosphere of Cologne, and took 
his doctorate in 1895 at the University of Bonn under the astronomer 
Friedrich Kiister. According to the Jahresverzeichnis . . . Universitéten 
. . . Abhandlungen, for 1894-95, his dissertation was printed under the 
name “Jean Peters.’”” Hence Deutsches Buiicherverzeichnis, Kiirschners 
Gelehrten-Kalender, Minerva, and other standard sources, preserved this 
form which was also copied on catalogue cards of the Library of Congress, 
and used by Klose, an associate of earlier years. From 1899 until his retire- 
ment in 1937 he was on the staff of the Astronomisches Recheninstitut of 
Berlin-Lichterfelde, where rigid formalism in name and otherwise was more 
the order of the day. Here he was the one principally responsible for the 
annual Berliner Astronomisches Jahrbuch, and Nautisches Jahrbuch. Besides 
elaborate tables for dealing with problems of pure astronomy (Kepler’s and 
coordinate transformations, for example) he published New Calculating 
Tables for multiplication and division by all numbers from one to four places 
(1909), an 8-place table (with Bauschinger) of the logarithms of numbers 
1-200 000 (1910), 7-, and 10-place (Anhang with J. Stein) tables of the 
logarithms of numbers 1-100 000 (1940, and 1922), Zweiundfiinfzigstellige 
Logarithmen (with J. Stein, 1919), the Factor Table giving the Complete De- 
composition of all numbers less than 100 000 (B.A.A.S., Math. Tables, v. 5, 
1935), and fifteen other volumes of tables of logarithmic and natural 
trigonometric functions. The first of these volumes was (with Bauschinger) 
the eight-place logarithmic table for each sexagesimal second (1911). This 
was later followed by tables of varying number of decimal places (3, 5, 6, 7, 
8, 10, and 21) for logarithmic, machine, and other calculations, with 
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sexagesimal division of the quadrant, with decimal division of the degree, 
with time as argument, and with centesimal division of the quadrant. Since 
during 1938-40 there were three editions of his six-place table for each 
0£.001, its possible use in connection with war problems is suggested. 

Klose refers to a large number of unpublished tables by Peters, but names 
only one, an eight-place table of reciprocals. Earlier reference has been made 
to eleven others; four of these (MTAC, p. 39, 114) being seven-, eight-, and 
ten-place tables of natural trigonometric functions for the centesimal 
division of the quadrant, and eight-place table of addition and subtraction 
logarithms (see L. J. Comrie, Br. Astr. Assoc., Handbook, 1929, p. 39). 
From information at hand it seems possible that one of these, the seven- 
place table for each 0*.001, was really published at Leipzig in 1939. The 
seven-place logarithmic table for each 0*.001 appeared in 1940; see RMT 125. 
Fifth and sixth unpublished mss. are referred to in RMT 78, 127, and five 
more in RMT 125. 

“‘J. Peters was an accurate and speedy calculator capable of sustained 
effort. He excelled in organizing great calculating tasks. His methods of 
calculation, and of avoidance of errors, were well thought through. Unfor- 
tunately these are scarcely known except to a small circle of his most 
intimate collaborators.” The excellent appearance of type-display in his 
published tables is also noted. 

[There are references to Peters and his works in MTAC, p. 7, 8, 10-15, 
25, 26, 39, 42, 43, 45, 55, 57-59, 64, 65, 82-84, 86-89, 100, 114, 121, 122, 125, 
143-148, 162-164, 170, 171]. 

R. C. A. 


14. RMT 98.—It seems worth while supplementing this review of 
MATHEMATICAL TABLES Project Miscellaneous Physiacl Tables, by re- 
marking that, in Table V if we regard x as a sine, G and xG are secants 
and tangents respectively. Similarly if x is a cosine, G and xG are cosecants 
and cotangents. This purely trigonometrical relationship is not pointed out 
in the Introduction. It is quite evident that this table may be used for 
applications other than electronic functions. x (review) = 8 (Table V). 


L. J 


15. Russian oF BaRLow’s TasBLes.—The third edition of 
Barlow's Tables of Squares, Cubes, Square Roots, Cube Roots, and Rectprocals, 
of all Integer Numbers up to 10000, edited by L. J. ComRIE, appeared at 
London, Spon in 1930 with xii, 208 p. In MTAC, p. 17, we noted incidentally 
the pirated Russian government edition of this. It has two title-pages; the 
one on the left intended to be the English original, while opposite it is the 
following (when transliterated): Tablitsy Barlou kkadratov, kubov, kornet 
kvadrainykh, kornet kubichnykh i obratnykh velichin iselykh chisel ot 1 do 
10 000. Gosudarstvennoe tekhniko-teoreticheskoe Izdatyeliestvo, Moscow 
and Leningrad, 1933. xii, 208 p. There were 10 000 copies in the edition. 
The last words of the title on the English-title page of the Russian edition 
are, “‘all Integer Numbers up 1 to 10 000.” 
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16. A Russian TABLE OF Propucts.—To meet the needs of the great 
army of workers on Collective Farms, in Industry, and in Factory Kitchens, 
and Cooperative Stores, in its second five-year period, the Russian Govern- 
ment published in 1936 an edition of 10 000 copies of the following bound 
volume: I. P. ZoLoTAREv, Universalnie Raschetnie Tabliisy. Prakticheskoe 
Posobie dlia Bukhgalterov i schetnykh rabotnikov. (Universal Calculating 
Tables. Practical Help for Bookkeepers and Computing Workers] Moscow 
and Leningrad, 565 p. +1 p. errata. 11 X 14.5 cm. Scattered through the 
volume are 15 other errata slips with one or two corrections on each. By 
means of this table one may read off at once the products of two numbers NV); 
and Nz, if N; = 0.01(0.01)0.07,0.10(0.05)1(1)250; and Nz = 1(1)250,255- 
(10)525. Examples of practical problems solved by the use of such a table 
are given on p. 3-9. 


QUERIES 


7. BricGcs’ ARITHMETICA LOGARITHMICA.—Compare N8. Since the publi- 
cation of J. W. L. Glaisher, Report of the Committee on Mathematical Tables, 
London, 1873, of J. Henderson, Bibliotheca Tabularum Mathematicarum, 
part 1, Logarithmic Tables, Cambridge, 1926, and of A. Fletcher’s letter in 
Nature, v. 148, 1941, p. 728 (to which my attention was drawn by L. J. 
Comrie), it has been generally known that there were some copies of Henry 
Briggs, Arithmetica Logarithmica, 1624, 388 p., with an extra 12 pages con- 
taining the logarithms of numbers 100,001(1) 101,000, to 14D, and the square 
roots of integers 1(1)200, to 11D, with first differences in each case. There 
are copies of such enlarged editions in the library of Trinity College, Cam- 
bridge, and in the Savile collection of the Bodleian Library, Oxford. Glaisher 
tells us that in 1873 C. W. Merrifield ! owned such a copy. Mr. Fletcher 
reported that the Harold Cohen Library of the University of Liverpool had 
a copy of just the 12 supplementary pages of the work of Briggs. Where are 
other copies of the Arithmetica Logarithmica, with the extra pages described 
above? 

R. C. A. 


1 The writer has many of Merrifield’s letters to Glaisher, 1873-81, mostly dealing with 
questions about tables. 


QUERIES—REPLIES 


7. AN ENGEL TABLE (Q6).—‘I am the owner of a copy of this table of 
natural trigonometric functions, containing 95 pages (19.7 X 29.9 cm.). The 
volume was an outcome of a meeting of German, Austrian, and Hungarian 
geodesists at Berlin in Nov. 1917, when the decimal division of the non- 
agesimal degree was agreed upon as desirable for tables. The Preussische 
Landesaufnahme undertook the preparation of the three volumes of 
logarithmic trigonometric tables necessary for this purpose, and published 
10-place (1919), 7-place (1921), 6-place (1921) tables, each for sine, tangent, 
cotangent, cosine, for every thousandth of a degree. These volumes! were 
all by the remarkable J. T. Peters. As a contribution from Austria to 
machine calculation ENGEL initiated the computation of the tables pub- 
lished in 1920, under the title given in the Query. It is a 10-place table for 
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sine, cosine and tangent, for every hundredth of a degree, with first dif- 
ferences. But additional tabular material makes possible the ready deter- 
mination of the values of the sine and cosine for each thousandth of a degree 
in the quadrant, and of the tangent for each thousandth of a degree 0° to 45°. 
Errors may occur in the values of sine and cosine up to 0.3 of a unit in the 
tenth decimal and up to 2.5 units for the tangent. For more exact computa- 
tions, a table of the interpolation coefficients for the second differences, 
progressing by 0.01, is added.’—This information was received by the 
EpiTors in a letter, dated 5 November, 1943, from L. W. PoLLaxk (66 Holly- 
brook Road, Clontarf, Dublin, Ireland), formerly professor of geophysics 
and director of the Meteorological Observatory at the German University 
in Prague. 

1 To meet the needs of an optical establishment Peters had earlier prepared a an 
table of these trigonometric functions for every thousandth of a degree; this was pub- 
lished in 1918 (see RMT 79). For a number of years thereafter, however, this volume was 


not readily accessible to the public. Then later (1935) came E. Buckingham’ s 8-place table 
for every hundredth of a degree (see MTE 12).—Epitors. 


CORRIGENDA 


P. 4, 8 (2), 52, 56, for TaLLQuist, read TALLOQvist. P. 22, 1. 4 from bottom, for 20, read 16.5. 
P. 36, 1. 14 for For, read ‘For. P. 39, 1. 8, for decimal, read decimals. P. 46, 1. 13 from 


bottom, read Sekamoto. P. 52, 1. 6 from bottom, read f (PT242—* 2 u). P. 53, 


1. 3, supply the factor |f®*)(£)| to the third member of (3). P. 56, . 10, read J. C. 
Adams. P. 69, 1. 12, for any Tables, read any other Tables. P. 70, no. 13, 1. 2, for 8D”, 
read eight places [should be 6D]”; 1. 5, for bei x, read bei’ x; for are the 41, not “35 
items”, read include the 36, not “35 entries”. P. 71, 1916, 23, 1. 2, substitute n for each x. 
P. 72, 31, 1. 3, for SD or 6D, read SS or 6S; 38, 1. 2, for A. N. Lowan, read R. v. MIsEs. 

P. 78, 1. 1, for 101 [J, L], read 101 [J, K, L]; for percentage of, read percentage points of. 
P. 78, 1. 30, read for Biological, Agricultural and. P. 85, 1. 5, read C. W. Merrifield; 
1. 13 from bottom, read J. de MENDIZABAL TAMBORREL. P. 86, in MTE 9, delete the 
“erratum” on p. 8, 1. 10, reported by GERTRUDE M. Cox. P. 89, first table, for 21, read 
25. P. 96, footnote 2, for Parfait, read Parfaits. P.97, N 8,1. 5, for 1560-61, read 1560/61. 
P. 108, 3A, for 1941, read 1940. P. 113, 1. 15, for 15D, read 16D. P. 121, 1. 15, for 17° 
39’ 50”, read 17° 39’ 40”. 
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EDITORIAL NOTICES 


The addresses of all contributors to each issue of MTAC are given in 
that issue, those of the Committee being on cover 2. The use of initials only 
indicates a member of the executive committee. 


Either the next number of MTAC, or the one following it, will prob- 
ably be a greatly enlarged special issue, devoted to a Guide to Tables in 
Bessel Functions. 


Volume I is to contain at least 8 numbers continuously paged. 
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